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ABSTRACT 


The perturbation method proposed by Professor Hui is 
described. The method gives exact solutions for the perturbed 
flow over both sides of a flat plate which is oscillating 
with small amplitude and frequency at large angles of attack 
in steady supersonic/hypersonic inviscid flow provided that 
the shock remains attached. Using the strip theory concepts 
these solutions are extended to study the dynamic stability 
in pitch of a flat, periodically oscillating wing or arbi- 
trary planform shape, at large angles of attack. Finally, 
Hui's perturbation method is extended to include the effects 


of upstream disturbances on a stationary wedge. 
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ie PRObuCcTILON 


This thesis deals with oscillating bodies in inviscid, 
steady or unsteady, supersonic/hypersonic flow. Its sub- 
ject matter may be divided into three main topics. 

The first topic is covered in Sections II and III.A.1,2 
and constitutes the background material for the other two. 
In Section II the Eulerian governing equations and boundary 
conditions for unsteady flow are formulated and the two- 
dimensional shock and expansion steady flow results are given. 
The basic elements of the linearized potential flow theory 
are also included in this section. In Sections III.A.1,2 a 
perturbation method proposed by Professor Hul is presented 
(Ref. 6,7,8].2 The method uses as a basis the assumption 
that the unsteady flow over an oScillating flat plate, with 
attached shock waves at an arbitrary angle of attack, is a 
small perturbation from the steady reference flow and, for 
small amplitudes of periodic oscillations, it gives closed 
form solutions for the flowfield quantities in the disturbed 
flow regions. 

The second topic is covered in Sections III.A.3,B,C 
[Ref. 5). In Section III.A.3 the closed form solutions over 


the upper and lower sides of the oscillating plate are combined 


1 professor of Applied Mathematics, University of Waterloo, 
Ontario, Canada. NPS, Department of Aeronautics, Visiting 
Professor in the period January-August 1980. 


hdl 





and the strip theory approximation is employed to calculate 

the necessary quantities in the disturbed flow regions over 

a flat wing of arbitrary planform shape. The in-pitch stability 

derivatives for the three-dimensional wing are then obtained 

in closed form and the results are compared with other exist- 

ing theories. In Section III.C a comparison of these results 

with linearized potential flow theory results is included 

while, in Section III.B, the fundamentals of the linearized 

theory, as applied to three-dimensional wings, are presented. 
Finally, the third topic 1s covered in Section IV. Based 

on the same perturbation method this topic introduces upstream 

unsteadiness in the flow and its effects on stationary bodies 

are obtained. The upstream unsteadiness is of a fairly general 

periodic form to give the solution in the flowfield generated 

by a formation of bodies provided that the body originating 

shocks are not crossing and the expansion fan regions are 

not overlapping. An extension to the case of oscillating 


bodies is readily possible. 
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HireewNolenoy INVEocclD FLOW THEORY 


To describe the fluid motion two methods are available: 
the Lagrangian method and the Eulerian method. In both 
methods the fluid is regarded as a continuum, l.e., its 
matter is assumed to be continuously distributed. 

In the first method the fluid is assumed to be divided 
into infinitesimally small regions called fluid elements or 
fluid particles. The so-called particle point of view is 
then adopted and description of the fate of each individual 
fluid particle is sought. To determine the unknowns asso- 
Clated with each fluid particle, e.g., its position coor- 
dinates, density etc., a system of equations 1s set up by 
applying to each fluid particle natural laws such as Newton's 
second law of motion and conservation of mass and energy. 
These equations are known as the Lagrangian equations of 
Plaid motion. 

Although the Lagrangian description appears to be a 
natural way to approach the problem of fluid motion, the 
Eulerian description is preferred in general since it gives 
more insight into the problem, it is much simpler and in 
most cases one is not interested in the fate of each individual 
fluid element but rather in the properties of the fluid ata 
certain point of the flow field at a certain time. 

In the Eulerian method attention is focused on the vari- 


ous points of the space filled by the flowing fluid and a 


ge) 





description of what is happening at each of these points, 
in terms of quantities of interest such as pressure, density 
and velocity is sought. The flow quantities of interest 
are assumed to be functions of space and time, i1.e., to be 
scalar or vector fields. Thus in the so-called field point 
of view adopted here the fluid flow is characterized by the 
fields of velocity, pressure, density and so on and a fluid 
element occupying a certain point at a certain time assumes 
for its properties the values that are appropriate to that 
point at that time. To solve for these fields a system of 
equations is again set up by using, as before, natural laws 
such as Newton's second law of motion and conservation of 
mass and energy. These equations are known as the Eulerian 
equations of fluid motion. 

Throughout this thesis the Eulerian approach is used. 

In the next subsection the major steps in deriving the 
Eulerian equations of fluid motion are indicated and the 
equations are presented in the form in which they. will be 


used later on. 


A. UNSTEADY EULER EQUATIONS 

The Eulerian equations of fluid motion, called hereafter 
Simply equations of motion, may be set up either in differen- 
tial form or in integral form. Furthermore they may be 
developed either from the point of view of a certain fluid 
region that contains the same fluid elements for all times 


(control mass approach) or from the point of view of a fixed 
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volume in space through which different fluid elements flow 
through (control volume approach). 

In what follows, the equations of motion are set up in 
the differential form from the point of view of an infinitesi- 
mal fluid region. 

In this derivation we will then naturally be involved with 
the calculation of the time rate of change of flow quantities 
as we follow the fluid element around, the so-called material 
derivatives of quantities. 

A physical interpretation of the material derivative and 
its components is briefly included. 

1. Material Derivative 

Consider a fixed coordinate system and a fluid ele- 
ment situated at point r at some time tt Let Q(r,t) denote 
some fluid property Q of interest (density, velocity, etc.) 
associated with the point r at time t. The fluid element 
Situated there (see Fig. 1.1) will assume for its corresponding 
property Q the value Q(r,t). In a short time interval At 
the element moves through a distance As = VAt where V is 
its velocity at r and t. The element will then assume for 
Q the value appropriate to its new position r+VAt at time 
t+At. If we denote this value as Q(r+VAt, t+At), the 


change of Q in the time interval At is 


AQ = Q(r+VAt,t+At) - Q(r,t) 


tearred quantities denote vector quantities. 


ales 





and the rate of change of Q following the element around, 
D( ) 





usually denoted by Dt is 
DQ _ jim QlerVet, t+At) - Q(r,t) 
Dt At 
At+0 


Expanding Q(r+VAt,t+At) in a Taylor series we get 


2 
O(rt+VAt,t+dt) = aiz,t) + (22) at + 28 ater e ia. 
= = 5) 
ret oe oe ae 
30 379 2 
+ (sp)_ Vat + (—5)_ (VAL) +... 
Pic eS r,t 


where s denotes distance in the direction of the velocity 
Wat point r and time t. 
Using the Taylor series expansion the rate of change of 


DQ 
©, Dt becomes 


ey (1.1) 


where the derivatives are evaluated at r and t. 

The total rate of change of any property Q is thus com- 
posed of two parts. 

To see the physical interpretation of each of these terms 
consider a flow field which is at any instant spacewise uni- 
form but varies from instant to instant and a second flow 
field which is steady but not uniform spacewise. Consider 
also a fluid element in these fields which in a small time 
interval 6t is moving from position r in the flow field to 


position r+Vét. 
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AS =VAt 






path 


(Gear yeeat (teat) 


feege L-l. Lilustration of local, convective and 
material derivatives. 


The change of a property Q of the fluid element moving 


in the first field described above is S25 and the rate at 


which the property Q is changing locally at the point r is 
=o This local rate of change is known as the local deriva- 
tive and is the first component of the total rate of change 
mieequation (1.1). 

The change of a property Q of the fluid element moving 
in the second flow field described above is ($2) vét. This 
change which-is called convective change is necessary since 
the element has to have an appropriate value for its property 
Q at its new position r+Vo6t. The rate of change of the 
property is (S2)¥v and is known as the convective derivative. 


It is the second component of the total rate of change in 


SPawacion (1.1). 


ik Ys 





The sum of the local derivative and the convective 
derivative as given by equation (1.1) is known as the total 
or substantial or material derivative. The last term is 
probably more descriptive since the derivative is constructed 
following a certain material element around. This term will 
be used hereafter. 

Recalling that <2 represents the derivative of Q with 
respect to distance in the direction of the velocity V we can 
write equation (1.1) in the following form 


— 
Dt 


| 
a} a> 
ao 


+ (e, -gradQ)Vv 


where a. is the unit vector in the direction of V or simply 


ae 
t 


+ V-gradQ ee) 


W5 
moO 
@ 


where Q can be a scalar or a vector quantity. 
If Q is vector quantity, i.e., Q = A the convective 


derivative V-gradA can be expanded using the formula 


MegqradA = s(grad(V-A) -¥ x curlA-Ax curl W 


- curl(V x A) +V(divA) - A> (divV) ] iS)) 


2. Equations of Motion 


We will now set up the basic equations that govern 
the unsteady motion of an inviscid, compressible fluid. We 
initially regard as unknowns the velocity field (Gee) Ene 


pressure field p(r,t) and the density field o(r,t). We 
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want to establish relationships between these fields by 
applying to a certain fluid element the basic laws of nature; 
Newton's second law of motion, law of conservation of mass, 
law of conservation of energy. 

a. Momentum Equation 

Let us consider an infinitesimally small fluid 
Miement situated at position r at time t (Fig. 1.2). If 
V and p are the velocity and density of the element at r and 
t and if édv denotes the volume of the element, then the mass 
and the momentum of the element are pdév and pdévV respectively. 
Let us also denote by F the total force acting on the element 
at time t. 

Newton's second law of motion which is applicable 
to any mechanical system states that "at any instant, the 
rate of change of momentum of a system is equal to the resultant 
of all forces that are acting on the system at that instant". 
By applying the above law to the fluid element considered 
and by noting that the rate of change of momentum of the ele- 
ment 1s simply the material derivative of the momentum we 


get 
D me vs 
De bP vv) = F 


Since the mass of the element odv remains constant the 


above equation becomes 


pév ae - F ieee 


iL) 








aguse 1.2. Fluid Element 


The total force F is the resultant of the so- 
called surface forces and body forces. 

The body forces are forces that act throughout the 
body of the fluid such as the gravity force. These forces 
will be assumed to be small and will be neglected in this 
thesis. 

The surface forces on the other hand are internal 
forces in the nature of actions and reactions across the sur- 
fact that separates the fluid element from its neighboring 
fluid elements. For a frictionless or inviscid fluid the 
Surface forces are simply pressure forces that act normal to 
the surface of the fluid element. Their resultant is 


- {Ff pnas 


6S 
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where p is the pressure of the element, 5S its total surface 
and n the outward unit normal to the surface. 

According to the integral definition of the gradi- 
ent of a scalar function the resultant of the pressure forces 


which is the total force F can be written 


a= -S[prnas = -6v gradp 
6 


Equation (1.3.1) then becomes 


— =  -gradp (1.4) 


oO 


This equation of motion is one of the fundamental equations 
of fluid dynamics and is called Euler's Equation. It repre- 
sents a system of three scalar equations for the five un- 
knowns--the pressure, the density and the three scalar com- 
ponents of velocity. 
b. Mass Equation 

Consider a fluid element situated at point r at 
time t with volume év (Fig. 1.3). The mass of the element 
1s opdv. The law of conservation of mass states that "the 
mass of any fluid element remains constant as it moves about" 
even though, in general, its shape, volume and density may 
change. 

Applying the law of conservation of mass to the 
element considered is thus equivalent to setting the rate 
of change of mass equal to zero or setting the material 


derivative of mass equal to zero. This gives 


ou 





or 
év = 0 (aro) 


The material derivative of the volume of the element may be 
expressed in terms of the velocity field. Let the surface 

of the element at time t be 6S and let it in a time inter- 

val é6t grow and become 6Sy as shown in (Fig. 1.3). The 

change in volume of the element 1s equal to the volume swept 

by the surface of the element during the time ét. If nis 

the outward normal to the original surface the net volume swept 


outward by 6S in time 6é6t is given by 


=] 
</ 


oe 
ids 
2S 


8S, 





Figure 1.3. Change in volume of a fluid element 


Zz 





GD Voenas 
6 


where V is the velocity. 
The material derivative of the volume of the 
element is then if we also employ the definition of the 


divergence of a vector 


Oe GV nas 2) fap Gidee a 
Dt fs 


Equation (1.5) then becomes 


O 


sr te div ¥ = 0 (1.6) 


This equation is known as the equation of conservation of 
mass, or simply, the equation of mass, or the equation of 
continuity. It is a relation between the velocity and den- 
sity fields only. Since it does not involve any dynamical 
quantities (such as pressures or forces) it is a kinematical 
relation. 
Gs Energy Equation 

The law of conservation of energy expresses the 
balance of energy exchanges that take place between a system 
and its surroundings. A fluid in motion may be regarded 
as a thermodynamic system characterized by the usual thermo- 
dynamic variables such as entropy, internal energy, etc. 
We will assume that the fluid is non-heat conducting and 
also that for a fluid element the only possible energy ex- 


change process is work done by the surface forces and body 
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forces. The law of conservation of energy applied to a 
fluid element may be expressed as follows. 

The rate of increase of energy E of a fluid 
element = the rate of work Wy done by the surface forces and 
the rate of work W. done by the body forces. Symbolically 


ee We +) W 


DE l (lives) 


Z 


We will assume, as before, that the body forces can be 
neglected and that the fluid is inviscid. The only possible 
energy exchange is then work done by the pressure forces. 


The rate of this work is 
W, = - (Yenas-V = - Of ev-nas 


where 6S as before, is the surface of the element, n is the 
outward unit normal and p, V are the pressure and velocity 
of the fluid element. 

According to the definition of the divergence of 


a vector we can write 


w, = -Qfpi-nads = -év div(pv) 
53 


On the other hand the energy E of the fluid element 
is the sum of its kinetic energy and internal energy. We 
specify the internal energy of the fluid by the scalar field 
e(r,t) which denotes the internal energy per unit mass at 


point r and time t. Then since the kinetic energy per unit 
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2 
mass is the total energy of the fluid element with mass 


odv is 


2 


Equation (1.7) can then be expressed as 


2 
(e +4) = = div(pv) cima) 


iS) 
S4[© 
ct 


Equation (1.8) is referred to as the equation of conserva- 
tion of energy or simply the energy equation. An alternative 


form of this equation is 
P ar = -pdivvV Gino) 


which can be found by subtracting from (1.8) the so-called 


equation of mechanical energy 


2 
(>) = - V-gradp 


te} 
=4|2 
ct 


The equation of mechanical energy is formed by multiplying 
both sides of Euler's equation (1.4) by V. 

The energy equation (1.8) has introduced the 
internal energy of the fluid as an additional unknown in the 
formulation of the governing equations. The list of unknowns 
includes the three scalar components of velocity, and also 
the pressure, density and internal energy of the fluid, while 


there are five equations available. 
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At this point we assume we are dealing with a 
perfect gas and introduce the equation of state for a per- 


fect gas 
Pp = op RT Gr. LO) 


aS a Sixth relation between the unknowns. 
We specify the temperature T of the fluid as a 
scalar field T(r,t) and since we are dealing with a perfect 


gas express the internal energy e of the fluid by the relation 


(=) 5 Came Clee 1) 


where C,, is the specific heat at constant volume of the gas. 
From equations (1.9), (1.10) and (1.11) we can 
find that the energy equation for a perfect gas may be 
written in the following form which will be used hereafter 
st‘) = 0 AeaT) 
p 
Summarizing we state that the basic equations that govern 
the unsteady motion of a non-heat conducting, inviscid, per- 
fect gas with constant specific heats are equations (1.4), 
(1.6) and (1.12). These equations are rewritten below for 


easy future reference. 


feontinuity) se 4 weou) =. O (la) 
Vs = We 7. 

(momentum ) net V-gradv + 5 0 (1b) 

(energy) Bee veyi-) = 6 (1c) 
dt 0 Y ot ) 
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fee) Boundary Conditions 

Physical conditions that should be satisfied on 
given boundaries of the fluid are known as boundary condi- 
tions. There are several types of boundaries and consequently 
there are various possibilities for the boundary conditions. 

We will consider two types of boundaries which are 
of main importance in this thesis, (1) "the solid-fluid 
boundary" where the fluid is bounded by a solid surface and 
(2) "the fluid-shock-fluid boundary" where two regions of the 
same fluid in different states of motion are separated by a 
flow discontinuity. The possibility of an infinitely weak 
discontinuity will not be excluded. 

The nature and number of the boundary conditions depend 
also on the form of the differential equations that govern 
the motion of the fluid. In this sense there are differences 
between the boundary conditions for a viscous fluid or an 
inviscid fluid. In the following the conditions for an invis- 
cid fluid are considered. 

a. Conditions at a Solid-fluid Boundary 

We assume that the fluid is bounded by an imper- 

meable solid wall and require that no fluid should cross the 
solid surface. Since the surface itself may be in motion we 
denote by V the velocity of the fluid and by Ve the velocity 
of the surface. The relative velocity between the fluid and 
the surface is ae Let the equation of the surface be 
given by S(r,t) = 0. A unit normal to the surface is then 


given by 
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———— rads 
Fioradsl 


and the component of the relative velocity normal to the 


surface is given by 


— 3 il = “=, 
eV A) n = E]gradFy <” gradS Mie grads ] 


Now assume that an observer moves with the surface particles 
that compose the solid surface. The observer cannot observe 
any change in the function S(r,t) considered as a scalar 
field. This means that the total rate of change of S(r,t), 


following a particle of the surface around, is zero, i.e., 


aS — | _ 
ae + Me gradS = 0 


The component of the relative velocity normal to the sur- 


face then becomes 


T-7)-n = tS 4 7. 
(V V5) n = FigradF] 3 + V-grads) (ey) 


The condition of impermeability of the solid surface is 
(V-Vo) cn = 0 
or 


36. ii _ = 4) _ 
DL = ara So AVS grads = 0 At Site te) 0 Cle) 


If the solid-fluid boundary is formed by the surface of a 


Stationary rigid solid the above equation reduces to 
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VegradS = 0 Reest)) = 0 (eS } 


It 1S pointed out here that the condition for- 
mulated above states that at each point of the solid-fluid 
boundary the normal to the surface component of the relative 
velocity between the fluid and the solid must vanish. Thus 
for an inviscid fluid nothing can be said about the tangen- 
tial to the surface component of the relative velocity which 
May or may not be zero. In short the so-called no-slip con- 
dition does not apply to an inviscid fluid. 

b. Conditions at the Fluid~-Shock-Fluid Boundary 

Consider one dimensional adiabatic constant-area 
Flow of a perfect gas through a discontinuity (Figure 1.4a). 
Assume that the flow quantities in regions 1 and 2 are con- 
Stant throughout the regions. The equations of continuity, 
momentum and energy between cross sections l and 2 give 


Meet. 1: pp. 55,56] 


“bea, lope, 
2 2 
Peco PT P54 
oe p #e p 


These equations hold as long as sections l and 2 are 
chosen outside the discontinuity region. The discontinuity 
region may be assumed to be vanishingly thin and sections l 


ema 2 may be brought arbitrarily close together. In this 
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Region | 


Bee? 


Figure 1.4a. Change of conditions across a normal 
shock in a constant area duct. 
Steady flow. 





Region Region 2 


Pagure 1.46. Change of conditions across an arbitrary 
ciscontinulity surface S(r,t). Unsteady flow 
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case the requirement of a constant area duct is dropped and 
the equations apply locally across the discontinuity. Con- 
sidering the discontinuity as a boundary separating regions 
l1 and 2 we may refer to equations above as boundary conditions 


at the discontinuity and write symbolically 


[pu] = QO (1.16) 
At x = xX 
2 O 
fopu' + p) = QO (7) 
ae ae = 0 (1.18) 
2 vy-lo 7 ; 


where the square brackets denote the change in the enclosed 
quantity across the discontinuity. 

A generalization of the simple one-dimensional 
flow problem considered leads to the well known Rankine- 
Hugoniot conditions in the form that will be used in this 
thesis. Consider two regions in space separated by a surface 
S(r,t). Assume that adiabatic flow of an inviscid, perfect 
gas is established from region 1 to region 2 (Figure 1.4b). 
Assume also that the surface S(r,t) represents a discontinuity. 
In general the equation of the discontinuity is not known 
a priori but will be found as part of the solution of the 
flow problem. Thus the boundary is a so-called free boundary. 

Let the flow quantities in region 1 be V, (x,t), 

P, (r,t), p, (r,t) and in region 2 Vi(r,t), p(r,t), p5(r,t). 
mae boundary conditions (1.16), (1.17) and (1.18) should apply 


locally at any point of the boundary S(r,t) provided that the 
\) 


Sb 





velocities u,, Us are replaced by the normal components of 


velocities ue ; Vie - For a boundary moving with velocity 

li 2 
Y the normal to the boundary components of the relative 
velocities between the fluid and the free boundary are given 


ey tOrmula (1.13) 


IL aS 


Ya = t1grads] ‘3e a ee) 
V = J (22 2. G -gradS) 
no +]grads | eae 2 


From equations (1.16), (1.17) and (1.18) using the same 


notation we get: 


At S = 0 
_ aS _ 
(continuity) [Lo apse ay Ve Sie ua 0 eles 9) 
(normal [o (= + Wow? = aie) oy eee L2G) 
momentum) 
1 (8S .F.y5)2 4K Piys)?] = 
(energy) [slap +V VS) ar BN ] 0 (file oe) 


where. the square brackets denote, as before, the change in 
the enclosed quantity across the discontinuity and the symbol 
"VY" stands for "gradient". 

The conservation of tangential momentum is not 
expressed by any of the equations above. To find the tan- 
gential momentum equation we require that the velocity com- 
ponent tangent to the discontinuity be continuous. The 


tangential velocity components are given by 


i 





<i 
il 

3] 
x 
< 


< 
il 
oy 
x 
< 


where n is the outward unit normal to the surface 


n= VS = {s_,S S_} in Cartesian coordinates. The condition 


Y 
of conservation of tangential momentum then becomes if the 


square bracket notation is used: 


(tangential we. Sivieee [wl 
momentum) S S S 
X y Z 


at S = 9Q 


where [u], [v], [w] denote the change of x, y and z components 
of velocity across the discontinuity. 

The equation above imposes two scalar conditions 
at the discontinuity, as one should expect from physical 
considerations. An alternative form of this equation is the 


following: 


(tangential 


erentum) [V xVS] = 0 At S = 0 (ive &» 


Memeeetons (1.19), (1.2)), (1.21) and (1.22) constitute the 
complete set of the Rankine-Hugoniot conditions. They are 
symmetrical and therefore remain unchanged if the brackets 
are taken to denote the change upstream rather than downstream 
through the discontinuity. A definite sense of flow direction 
is provided by the second law of thermodynamics, which requires 
that the entropy shall not decrease across a discontinuity. 


The change of entropy is given by [Ref. 1: p. 60]: 
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and the requirement stated gives the following condition 
which should accompany the Rankine-Hugoniot conditions 


(2nd law of pe 0 a3 2G = (1223) 


thermo.) al _ 


It should be noted that so far we have avoided using the 
term shock instead of discontinuity, though shocks are the 
only possible physical discontinuities. This was done on 
purpose since we intend to use the Rankine-Hugoniot condi- 
tions across hypothetical expansion fronts (negative or expan- 
Sion shocks) through an iterative procedure so that in the 
limit condition (1.23) should not be violated. 
4. Shock-Expansion Flows 

Unsteady flows with shock waves or expansion waves 
are considered in this thesis. These unsteady flows will be 
solved by first finding the corresponding steady flow solution 
and then using it as a reference flow in calculating the un- 
steady perturbation flow. 

Since a number of exact steady flow conditions are 
Mmieeady available, they will be utilized as reference flows 
in finding the governing equations and boundary conditions 
of the corresponding unsteady flows. They include the super- 
sonic uniform wedge flow and the Prandtl-Meyer expansion flow. 
The results for these steady flows are stated below for easy 


reference. 
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a. Wedge Flow 
Consider steady uniform supersonic flow past a 
symmetrical two-dimensional wedge with semi-vertex angle 9 
(Figure 1.5). The wedge is assumed stationary. Oblique 
shocks will be formed at angles 8 measured from the free 
Stream direction. By conservation of momentum the tangential 
component of velocity is continuous across the shock so that 


V =V, . Then V and V are related by the normal shock 
i «2 Dy sy, 


relations. Since V = V, sing and V = V. sind the normal 
ny L no Z 


Shock relations can be used directly with M, replaced by 


IL 


M, sing and M Sing. The resulting relations 


1 Z 
for the oblique shocks are 


replaced by M. 


P 2M 


9 l Sanck = eb 


P y+1 





Figure 1.5. Steady supersonic flow past a stationary 
wedge 
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D 3 (y-1)M%sin“B +2 
- 2yM%sin“8 = §O~ ah) 





Oe 2 
PF . (y+1)M,sin 8B 
Pi 2 + (y-1)Mfsin“s 
—2 _ tan B 
Py tan 


p 
From the last two relations eliminating = and recalling 
1 


that ¢ = 8-8 we get 


Misin“s “il 
tan 8 = GOums YET ie? — (w2sin“@ - 1) Cis 2 4) 
(My sin 6 -1L) 


One way to solve this equation, i.e., to find 8 for given 

M and 8 is to express it as a cubic equation in x = cot 8 and 
select the appropriate positive root corresponding physically 

to the weak shock wave [Ref. 2: pp. 452-453]. The following 
equivalent equation was used in this thesis for numerical 


calculations 


a+Yely SPeNeeycone x? +(1 +421 wv? 


2 _ 


where x = tan 8B. 
For attached shock waves this equation gives three 
real roots for 8 and the middle one is the one corresponding 


to the weak shock wave. 
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A maximum deflection angle? for given M, (beyond 


1 
which the shock becomes detached ) can be found from equation 
mee-) by differentiating it with respect to 8 and equating 
eo) Zero. 
b. Prandtl-Meyer Flow 

Consider steady, two dimensional, wniform super- 
sonic flow over a convex corner (Figure l.6a). A turn of 
the flow through a single oblique expansion wave is not 
possible since this would lead to a decrease in entropy. 
The flow expands isentropically through an infinite number 
of centered straight Mach lines that form the so-called 
Pranctl-Meyer expansion fan. Thus upstream of the ray OB 
where @ = @. the flow is uniform with Mach number M and 
downstream of the ray OC where 6 = o4 the flow is also uni- 
fOrm with Mach number M,. For angles 6 such that @. < @ < 8 


a 
the flow field has the same properties along any ray 


L 


8@ = constant. The polar coordinate system shown has been 


TT 


chosen so that [lf =a + > + Biewnere er (Mj eis the Prandtl— 
Meyer function 

P(M) = = wa i eer - vem wean eee 
with A = fey +72 and y the ratio of specific heats of the 


y+1 


gas. The Mach number M, 1s given by 


tor the wedge at an angle of attack a the flow deflec- 
tion angles are @+a and 6-a for the lower and upper sides 
respectively. 
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Magure L.ba. Steady Supersonic flow over a convex corner 





“=stagnant fluid region 


Figure 1.6b. Expansion over a corner with a > I(M) wax 
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Along any ray 6 in the expansion fan the Mach number pressure 


and density are given by the following relations. 








@ = = tan™* [a (m*-1) 777] 
E(M,) 7 
E(M_) 
2 o° »1/ (y-1) 
where: 
E(M) = 1+ Y= ue 


The r- and 6-velocity components in the expansion fan are 


given by 


V = Cc sin 6 


< 
II 


Ac cos 86 


where: 


It can be seen that the Prandtl-Meyer function 
P(M) defined above has a maximum value for M> ~. This value 


7S 


S\2) 





— ui “ io = 
os) ee a 5 (V (y+ DAy 1) 


For convex corner angles a > P{M) 5 - P(M_) the streamlines behind 


x 
the expansion fan behave as though the flow occurred over an 
expansion of EA) ee and for an inviscid fluid a region of 
stagnant fluid lies between the hypothetical position of the 
body as sensed by the flow and the actual position of the 
body (Figure 1.6b). Whenever this occurred in numerical 


calculations performed, the influence of the flow on the sur- 


face of the corner! was assumed negligible. 


B. LINEARIZED POTENTIAL EQUATION 

The so-called linearized theory of supersonic flow builds 
up the flow produced by the motion of a body by superposi-~ 
tion of small disturbances such as those produced by a moving 
sound source. One can develop in this way relatively simple 
methods for the computation of velocity and pressure distri- 
butions in the field. 

In the case of vortex-free flow, the equations of motion 
can be reduced to equations analogous to the wave equation. 
The coordinate parallel to the direction of the main flow 
plays the role of the time coordinate. Hence the methods 


of finding solutions of the wave equation can be used. 


lie mean that in the equivalent case of a flat plate 
at an angle of attack a the pressure on its upper surface 
was assumed Zero. 
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The linearized theory however has serious limitations. 
First, it gives only a first approximation since all devia- 
tions from the uniform parallel flow are considered small 
and therefore additive. This is justified for very thin or 
slender bodies at small angles of attack only. Second, there 
are speed ranges in which the linearization of the equation 
of motion even for small disturbances is not justified. For 
the linearized theory to be valid the following two conditions 
must be met. 

(a) The perturbation velocities must be small in compari- 
son to both the main stream velocity and the velocity of sound. 
This condition excludes the case of very high velocities 
since if the mean stream velocity is several times larger 
than the sound velocity, disturbances which are small rela- 
tive to the Rana Steeen velocity may be of the same order of 
Magnitude as the sound velocity. This speed range is called 
the hypersonic range. 

(b) The perturbation velocities must be small in comparison 
to the difference of the main stream velocity and the sound 
velocity. This condition excludes therange near M = 1, the 
so-called transonic range. 

In spite of the limitations described above the lineari-~ 
zation of the equations of inviscid, compressible fluids 
proved to be of excellent use in developing approximate solu- 
tions in the supersonic range. 

There are three general methods used in the linearized 


theory of supersonic aerodynamics. 


41 





(a) The method of fundamental solutions or sources. This 
method is based on the superposition of fundamental solutions 
of the linearized hyperbolic equation for small perturba- 
tions of a uniform ae flow. In formulating wing 
problems this method uses sources and doublets located in 
the plane of the wing and the strength of these singularities 
is determined so that the boundary conditions applicable to 
the wing planform and shape are satisfied. 

(b) The methods of acoustic analogy and operational cal- 
culus. In these methods the solution of the hyperbolic equa- 
tion is expressed by means of Fourier and Laplace integrals 
respectively. The second method is better adapted to super- 
sonic flow problems since the Laplace integrals exclude the 
possibility of upstream travelling signals while in the case 
of Fourier integrals one has to impose additional conditions 
to secure that this possibility is excluded. 

(c) The method of conical flows. This method is based on 
conical flows, i.e., flows for which the velocity components 
at points lying on a straight line drawn from a point chosen 
as vertex are independent of the distance from the vertex. 

In this method the solution of a hyperbolic equation in three 
variables is reduced to the solution of Laplace's equation or 
wave equation in two variables and the existing methods of 
conformal transformations and the theory of. functions of com- 
plex variables can be employed. 

Methods of higher approximations, i.e., methods which 


lead from the simple case of the linearized solution toward 


\) 
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the exact one in successive steps, thus extending the range 
of satisfactory approximation, are required whenever the 
perturbation velocities are not small compared to the main 
stream velocity. 

In this section the concept of irrotationality will be 
introduced and then the basic assumptions and steps followed 
in deriving the linearized potential equation and the appli- 
cable boundary conditions for the general flow problem of a 
body in supersonic flow will be described. In the next sec- 
tion the linearized wing problem will be described. 

ie Lrrotational Flows 

By potential flow we mean that the velocity V is 
derivable from a scalar velocity potential 9, i.e., v = grado. 
On the other hand the vorticity or rotation w in a fluid is 
defined as w = curl V and the flow is called irrotational 
mee = 0 or equivalently if 

OV du OW OV du dW 
ax Jy Jy OZ dZ OX 
where u, v, w are the x-, y- and z- components of V. 

Physically the irrotationality of the fluid means 
that the fluid particles have zero moment of momentum about 
their own center-of-gravity axes or simply that they remain 
parallel to themselves as they move around. 

The condition of irrotationality of the flow is a 
necessary and sufficient one for the assumption of a poten- 


tial flow since the mathematical identities curl grado = 0 
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and div curl V = 0 show that the velocity V can then be 
Smmeewin the form V = grado. 

The irrotationality throughout a flow field can be 
proved by applying the theorems of Kelvin and Stokes. 
Stokes' theorem states that the area sum of the rotation over 
a given area is equal to the integral of the velocity around 


a curve bounding the area. Formally, 
jeerdiw— (curl V-dA = DV-dy 


The line integral is called circulation and is denoted by 
rT. Thus Tf =Q V-d& = fw-daA. 

As a consequence of Stokes' theorem, w = 0 if the 
circulation [ vanishes for all paths wholly within a simply 
connected flow region. 

Kelvin's theorem on the other hand states that the 
circulation [ about any contour always composed of the same 
fluid particles (i.e., a fluid line) is constant in an inviscid 
fluid with only conservative or irrotational body forces. 

ONE 


Mer an inviscid fluid it states that oa -$2 and it reduces 
pl 


EO De = QO when there is a simple relation connecting p and 
op. Physically it means that circulation [ about any line 
contour remains constant in time as we move along with the 
mruid. 

It is a consequence of both theorems that initially 


irrotational flows originating in a reservoir under uniform 


stagnation conditions or from straight parallel streamlines 


as 





will remain irrotational throughout the flow field at all 


times if there are no shock waves. 
Pmeeeinearized Potential Flow Equation 
In view of the above reasoning we assume that the whole 
flow field is irrotational and set V = grad®#. As a result 
the three unknown velocity components u, v, w are expressed 


in terms of the scalar field ¢% by 


By using relation (1.3) the momentum equation (1.4) 


may be written as 


oa 2 
0 (25 + grad = - Vxcurl V) = -gradp 


or since curl V = 0 and v = grad we may write 


Z 
3 ® V dp = 
grad (~— ar ay + f ) = 0 


By integration we get 


FG) 


Q?} a 
rt, o 
+ 

| < 
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If we define 6 = ®-fF(t)dt this equation becomes 


yA 
d ts 
Pb Spe fee = 0 


Q 
& 


Q 
ct 


By differentiating this relation with respect to time 


and by taking its gradient we get the following two relations, 
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if the additional relation 0,2 = Se for the velocity of 
sound is used. 
2 2D 
moog we) 6. 6 OR 
ot ot 2 On out 
2 2 
9 ¢ V _ a 
- grad (a os =) ice grado 


We now introduce continuity equation (1.6) which, since 


ay VY = 779, may be written as 


00 a 


5 grado + 7~6 = 0 


Ol[F 
ols 


Introducing in this equation the last two relations from 


momentum equation we get 


26 yyv2 
dt 


Z 1 .da 
yo - S{25 + 


pia ye 
5 + V grad (->) | = Q Gls 5:} 
Gy ot 


Equation (1.25) is the exact non-linear differential equation 
to be satisfied by the velocity potential % for an unsteady, 
inviscid, irrotational flow. Because of its strong non- 
linearity, solutions have been found in very few special 
cases. Thus the small disturbance concept is introduced which 
leads to linearization of the equation. 

We assume that the velocity vector V differs only 
Slightly in direction and magnitude from the free stream 
velocity U_, taken along the x-axis, and we define a distur- 
bance velocity potential ¥ obtained from the total velocity 
potential % by separating out the contribution of the uniform 
flow. In this way we have 
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The local velocity components are then given by 


t— U + u', v = v', w = w' 


where u', v', w' are the perturbation or disturbance veloci- 


ties which can be found from 


We regard all disturbance velocities small in compari- 
son with U,, a and U_-a and all pressure and density changes 
small in comparison with main stream pressure and density. 

We also assume that the small quantities change gradually 
in all directions and that time variations are not too rapid. 

We now return to equation (1.25) and assume that 
the linear terms in Laplace's operator are of the same order 
of magnitude. Substituting the velocity V in this equation 


and retaining first order terms only we get 


2 2 Z 
2 2 hae ae oY Zo ¥ - 
y =) aes + 2U., sae UL 5! = 8 (1.26) 
Cuero & ox 


where the second term in parenthesis is found from 


] 


IL = _ di: t = =- i 
Meese) = = =sil2(u,tu')¥.1 = pete te 


O Oo io 


and the third term in parenthesis is found from 
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To complete the linearization of (1.26) we should dispose 


of the factor —. 
a 


For steady flow ie may be eliminated by using the 
following relation which 1s a consequence of the first law 


of thermodynamics 
a +—-= V = constant 


For unsteady flow let a slender body move in the flow 
field and denote the velocity of sound far upstream in the 
Maaesturbed field by the constant a. In the vicinity of the 

2 


body a will then be a variable depending on position and 


can be represented by a sum of terms of the form 


at = a’ + (Aa) * + ... 


As a first approximation we set a = a, and introduce this 


value in equation (1.26) getting 


2 2 2 
2 lL .a“Y¥ 9“ 2 
ee 2 = gxse 9. 2 = aa 
oe cae dX 


CG 


An alternative form of this equation with M, = — is 


co 


2M 


2 = 1 - 
(1 Me) ax t yy F422 7 OD Fs Fee = 0 (1.28) 


Equation (1.28) is the linearized unsteady potential equa- 


tion and is used as the basic equation in most aerodynamic 
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analyses. Because of the assumptions made and terms retained 
this equation is valid for unsteady, inviscid, irrotational 
Flows that are purely subsonic or purely supersonic and is 
limited to small disturbances only. 
3. Linearized Boundary Conditions 

To completely specify the mathematical problem that 
describes the flow, the following boundary conditions need 
generally be prescribed [Ref. 3: pp. 1.27-1.29]. 

a) Surface Boundary Condition: The wing surface is 
impenetrable to the medium. 

b) Edge Conditions: Enough viscosity remains in the 
inviscid fluid to determine the flow pattern near sharp 
edges. 

c) Wake Conditions: The free vorticity shed from the 
trailing edge must have a circulation which vanishes together 
with the bound circulation. It is furthermore assumed that 
the shed wake is a continuous sheet of discontinuity which 
1s coplanar with the wing projection in the direction of 
flight. Edge effects and rolling up of the sheet are 
disregarded. 

ad) Conditions at Infinity: A state of uniform flow must 
be prescribed at inifinity. In addition the Sommerfeld 
radiation condition requires waves to propagate away from 
sources of disturbance toward infinity. 

e) Other Conditions: As the most important additional 


condition, the requirement that proper account be taken of 


49 





zones of influence and action at supersonic flow velocities, 
is mentioned. 
The first of the above conditions, namely, the con- 
dition of impermeability, as applied to a wing or airfoil, 
is considered next. The procedure can be easily extended to 
cover other bodies of interest such as slender fuselages, etc. 
Consider a wing fixed relative to a Cartesian co- 
ordinate system so that it lies close to the xy-plane (Figure 
1.7). Assume that the wing is submerged in an infinite mass 
fameeeid MOVing with velocity U, in the positive x-axis. 
Let the upper and lower surfaces of the wing be expressed by 


equations 


il 
© 


S = Zo Z (xy ,t) 


il 
ro) 


S = Z@- Z, (X,y,t) 


The condition of impermeability of these surfaces, 


by recalling equation (1.14), requires that 


dZ dz OZ 


= u u _u — } 
—- 3. + u a5 + Vv ay Ione vA Ze (Xen Ro 
92, 92 92 
— eS es ——— = } R 
Ww aoe aa tC e Oren Z Zor Cy) meen R, 


where u, v, w are the components of velocity V and R, is 
the portion of the xy-plane covered by the projection of the 
planform. 

These are exact, non-linear equations. To linearize 


them assume, as before, that the disturbance velocities are 
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mathematical plone 
surface 


Figure 1.7. Linearized boundary conditions for a wing 


small compared to the free stream velocity U_ and also that 


dz dz 
the slopes — aa etc., are very small compared to unity. 


Then retaining first order terms only we get 


dZ OZ 


= | a oa = 

W yt a ee LOL Z Za! (ae, Ve) Th Ro 
IZ OZ) | 

W = Syre or UL Ss fe te 2 = Si) G (oa, 7) E39) Ro 


Since Zz and z, are small compared to the wing chord 


Me 
we may, as a further step, replace the actual wing with an 
infinitesimally thick surface of discontinuities in u, v, w 
and pressure p. With this mathematical plane surface located 


on the xy-plane, we may expand w in Maclaurin series about 


its values just above and below the xy-plane 


oye 





> 
mee, Y , Z re) = w(x,y,0',t) ae ow (xX, ae) A) a 


u OZ 


W(X,Y,2,,t) = w(x,y,0 ,t) +z Eee ae ease “ 


Using the same arguments as before the higher order 
product terms can be neglected and the impermeability or 


flow tangency conditions take the following linear forms 


+ OZ, OZ 

w(x,y,0 ,t) = ar or Ee ce cay) in R. (1.29a) 
az IZo OZ) 

max, Ys 0 Bit) = =) = ly “Ia ce 7) in A (1.29b) 
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III. OSCILLATING WINGS OF GENERAL PLANFORM 
IN SUPERSONIC/HYPERSONIC FLOW  +~— 

Consider waco. (Sspacewise and timewise) supersonic 
or hypersonic flow of an inviscid perfect gas with constant 
Specific heats past a flat wing of a general planform shape 
at an angle of attack a. Assume that the wing is performing 
a small amplitude slow pitching oscillation. 

The problem considered is to find the unsteady flow 
quantities in disturbed regions over the wing and thus its 
stiffness and damping derivatives. 

The governing equations of motion are given by equations 


(la-lc), restated below. 


sc 7 V- (ov) = 0 (2.1a) 
yey ae = 0 (2.1b) 
at O 

ea eee ye) = «CO (Zc) 
at oY su 


where p, 9, V and y are the pressure, density, velocity 
and adiabatic exponent of the gas. 
The flow tangency condition to be satisfied at the sur- 


face of the body is given by equation (1.14), restated below. 


V-vVS = 0 At Ss = 0 (G22) 


= 
+ 
< 


where S(r,t) = 0 is the equation of the body surface. 
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The fluid~-shock~-fluid boundary conditions to be satis- 
fied across the shock are given by equations (1.19-1.22), 


restated below. 


[p(s + F-vG)] = 0 At G = 0 ey 
i a6 murvc)- (9G) -] = 0 (2.3b) 
[s(22 + F-va)* + 1, Bive)*] = 0 (2.3c) 
he<vVG} = O (225d) 
where G(r,t) = 0 is the equation of the unknown shock shape 


and the square brackets denote the change in the enclosed 
quantities across the shock. 

Equations (2.1-2.3) are nonlinear. The nonlinearity of 
the governing equations and boundary conditions along with 
the existence of a shock with an unknown shape, contribute 
to the complexity of the problem considered. 

For low supersonic Mach numbers and very low angles of 
attack shock waves can be replaced by Mach waves and the 
linearized supersonic potential flow theory can be employed. 
The problem can then be solved, at least for certain groups 
of planform shapes, by fairly general methods. The fundamen- 
tals of the linearized supersonic potential flow theory as 
applied to a three-dimensional oscillating wing are presented 
maesection IIi.B. 

For high angles of attack and/or Mach numbers the shock 


waves become strong and the linearized theory cannot be used. 
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To overcome the Aes euleies encountered in this case, namely 
the existence of the shock and the nonlinearity of the equa- 
tions, Professor Hui proposed the use of a perturbation method 
in which the unsteady flow field is regarded as a small 
perturbation to some reference steady flow. Thus the unsteady 
flow problem is to be solved by first finding the corresponding 
steady flow solution and then using it as a reference flow 
in calculating the unsteady perturbation flow. The solution 
of the three-dimensional wing by this method will be presented 
maesection III.A. 

Finally, in Section III.C results for the stability deriva- 
tives are presented. A comparison with linearized potential 


flow theory results is also included. 


Be PROFESSOR HUI'S THEORY 

In this section the problem of dynamic stability of a 
flat wing of a general planform shape at arbitrary angles of 
attack in steady supersonic/hypersonic flow is considered. 
The wing is assumed to be oscillating in pitch with small 
amplitude and frequency and the bow shock be attached to the 
body at all times. 

The problem is covered in [Ref. 5] and only the basic 
steps will be included here, in Section III.A.3. Its solu- 
tion is based on the assumption of an inviscid perfect gas 
with constant specific heats and the perturbation method 
developed by Professor Hui is employed to calculate the 


resulting unsteady flows over the upper and lower surfaces 


5 





Of a two-dimensional flat plate. Finally, the strip theory 
approximation is utilized to combine the effects of these 
flows for the case of a three-dimensional wing, provided that 
the bow shock is attached and therefore the flows are independent. 

The unsteady flows over the lower and upper sides of a 
two-dimensional flat plate are studied in [Ref. 6,7] and 
[Ref. 8] respectively. Nevertheless, we will indicate, in 
the following first two sections, III.A.1 and III.A.2, the 
way in which these flow problems are formulated and solved. 
We will also give the solutions for the complete set of flow 
quantities in the disturbed regions. These flow quantities 
will be used in Section IV where the effects of upstream 


unsteadiness in the flow are considered. 


1. Two-Dimensional Oscillating Flat plate -~-Compression 


Side 

Instead of a flat plate the equivalent flow problem 
of a two-dimensional wedge is considered. This problem is 
formulated and solved in (Ref. 6,7] with the ultimate goal 
of studying the stability of wedges/caret wings. In what 
follows in this and the next subsection, 

a) The major steps in the method of solution are 
indicated. 

b) A generalized approach that permits the formulation of 
the fluid-shock-fluid boundary conditions is adopted. This 
approach is described in Appendix A and the formulation of 
the boundary conditions for the two cases is given in Appen- 


dices B and C. 
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c) Some of the results and discussions contained in 
[Ref. 6,/7,8], pertaining to the stability of wedges/caret 
wings and flat plates, are not included since they are not 
directly related to the subject matter of this thesis. 

d) On the other hand the solutions for the complete set 
of the flow quantities in the disturbed regions, which are 
not included in the above references, are given. Much of 
the mathematical detail in obtaining these solutions is 
omitted. These flow quantities will be used, as mentioned 
before, in Section IV where the effects of upstream unsteadi- 
ness in the flow are considered. 

e) The same symbols as those used in the references will, 
in general, be employed. Changes will be limited. to those 
necessary for clarification purposes or generalization of 
approach. 

a. Problem Formulation 

Consider a two-dimensional wedge of length 2%, at 
design condition (zero mean angle of attack), in a supersonic/ 
hypersonic, uniform, steady flow of an inviscid perfect gas 
with constant specific heats (Figure 2.la). Assume that the 
wedge is performing a low amplitude and frequency harmonic 
oscillation in pitch with given circular frequency w, about 
an axis perpendicular to the plane of the paper, through the 
point C shown. Let a system of cartesian coordinates Oxy 
be attached to the wedge so that O is at its apex and axis 
Ox is along the mean position of the upper surface. The bow 


shock is assumed to be attached to the body and the flow 


5) 7 
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Figure 2.1b. Oscillating flat plate--Compression side 





quantities on the upper surface of the wedge are to be 
found. For all practical purposes the solution to this 
problem will give the flow quantities on the lower side of a 
two-dimensional flat plate of length 2 sec 9 at an angle of 
@erack 9 (Figure 2.1b). 
b. Method of Solution 

The unsteady flow over the upper surface of the 
wedge will be found by perturbing the steady shock flow con- 
Sidered in Section II.A.4.a. 


Denote by U., p., 9p, the velocity, pressure and 


oO 


density in region A. Denote by Ur Por P the velocity, 


O 
pressure and density of the steady reference flow in region 


B. Non-dimensional lengths and time are introduced, defined 


by 


| 


and. t= (2.4) 


* 

I 
| | 

a 

I 
=a 
=I Jo” 


Assume that, as a result of the oscillation of 


the wedge, the perturbed flow quantities in region B are 


given by 
oS One Bt gg (22 5.) 
Vos evt+.. (2. 5b) 
oa P, + Ep + ... (56) 
o = Q, + a6 fers i235) 
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a 


where « describes the deviation of the unsteady flow from 
the reference steady flow. 

Substitute quantities (2.5) in the governing 
equations of motion (2.1), non-dimensionalize the indepen- 
dent variables, by using relations (2.4) and get the following 


perturbation equations 








il 
u, +u = - p (2.6a) 
te Oe P44, x 
es + eee: =. p (2 0b) 
aa 
P, + Pp = a (6 Oc} (256) 
fe pid Or re x i 
Po 
et PL 7 ux a = 0 (2.6d) 


where subscripts denote partial differentiation and Oo is the 
speed of sound in the reference steady flow. Assume that 


the perturbation quantities u, v, p and op have the form 


u = ue **u(x,y) (229 a) 
7 ikt, 

v= ue (x,y) (2 4975) 
7 ete 

Deel P (x,y) erie 
a ec 

ee OMe BUC Gy) OR Esh) 


where U, V, P and R are unknown quantities to be found, Mo 


is the Mach number in the reference steady flow and k is 
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mae SO—-Called reduced frequency of oscillation defined by 
k = w/a. k is assumed to be small and the time indepen- 
dent quantities U, V, P and R are expressed as power series 


met ik) of the form 


(0) (1) 


U = JU eae) U + . (2.8a) 
may 8! ee (ikj)v'*? + (2.8b) 
meee 8) Gkj)p) +. (2.8c) 
peek? e GikyR’*’ + .. (2.84) 


Expressions (2.8) are substituted in (2.7) and the resulting 
expressions in the perturbation equations (2.6). By equating 
the same order terms, in each of these equations, a sequence 
of systems of partial differential equations is formed. 

Only the systems of zeroth and first-order equations are of 
interest in stability analysis. The zeroth-order equations 


are 


(Oo) _ 3 (0) 

u, = - = Pl (2.9a) 
O 

(Oo) 5 EG | 

Vy = i By (2.9b) 
O 

(a) eee) (2.9c) 

~ x 

y§ 4 yl 2, wy pl) 2 9 (2.9a) 

x y omex 
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The first-order equations are 


‘ " wo Py = -y'°) (2.10a) 
O 
y (4) oe p (1) _ _y (0) (2.10b) 
Xx Ms Y 
p (1) 7 pit) _ R69) = p (0) (2.10c) 
ut) yg vy) gw RO) 2 Ly Rf) (2.104) 
x y Omx O 


Next we consider the boundary conditions applicable 
to the problem. 
Along the surface, the condition to be satisfied 
is the flow tangency condition given by equation (2.2). 
The equation of the surface for a Sion nee wedge is given 
by S = y = 0 and for an oscillating wedge is given by 


y te (h cos 8% -x) = 0 where € = age Equation 


max, y , t) 


eee then, gives with V = {u,v} = uti +eU,evi, 
V(ix,y) = 1 + (ik)(*h cos 6) at y = e€(x-h cos @) 


Expanding V(x,y) about y = 0 and neglecting higher-order 


terms we get the linearized condition 
V(x,0) = 1+ (ik)K-hcos?) at y = 0 
Use of equation (2.8b) gives 


y 6) = 1 At y = 0 (2.9e) 
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vit) = x-hcos@ at y = 0 (2.10e) 
Across the shock, the conditions to be satisfied, are the 
conditions given by (2.3). The equation of the shock in 
the steady reference flow 1S given by Go = -y + x tan > = QO. 
Let the equation of the shock in the unsteady flow be given 
by G = -y + x tan > + €Q(x) = 0, where Q(x) 1s an unknown 
function to be determined as part of the solution and 
Ee = eer kt. To find the boundary conditions across the shock 
we substitute expressions (2.7) in (2.5) and the resulting 


expressions in equations (2.3). The boundary conditions, 


after linearization, are given by 


We AQ' + (ik) BO AT y = x tan 6 (2einica) 
P = cq! + (ik) DO (27.1) 
U = EQ! + (ik) FQ (25 bie) 
R = GQ' + (ik) JQ (2.114) 
where Q' = dQ/dx and the values for the constants A through 


~~ 


J, which depend on the reference steady flow, are given in 
Appendix B. The derivation of these relations is lengthy 

and tedious, even for the case considered here, where there 
are no upstream disturbances. In Appendix A equations (2.3) 
are put in an alternate form. This form permits a much easier 
solution of the equations and is repeatedly used throughout 
this thesis. " 


63 





Expressing the time independent quantities, U, 
V, P and R by their expressions as given by equations (2.8) 
we get the following boundary conditions for the zeroth and 


first-order systems respectively. 


At y = x tan y $9) - ag! (9) aot) 
Eo e=tecoll (2.9g) 
GO a CO (2. 9h) 
R'9)  . Go: (9) (2.9i) 

at y = x tan $9 vit) = agt (1) + Bq 9) (Zeon) 
P|) 2 egh |) 4 59!) (2.10g) 
Ul) Logg (1) 4 59 (0) (2.10h) 
pit? = got 1) 4 gQ (9) (2.104) 


c. Complete Solution 

Two boundary value problems have been set up. 

The zeroth-order equations (2.9a-2.9d) and the 
zeroth-order boundary conditions (2.9e-2.91) constitute the 
first boundary value problem. This problem, which will be 
solved first, corresponds to the . problem of steady flow past 
a wedge and its solution should give the flow: quantities | 
behind the shock for a stationary wedge. This result is 


shown in Appendix B. 
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The first-order equations (2.10a-2.10d) and the 
first-order boundary conditions (2.10e-2.101) constitute 
the second boundary value problem. 

In bothproblems the equations and boundary con- 
ditions are linear and therefore suggest for the unknowns 
i, as p+) and oe (1 = 0,1), solutions that are 
linear combinations of the non-dimensional spatial coordinates 
x and y. 

In view of the above reasoning, we assume for the 


first problem a solution of the form 


0 0 0 0 

yl _ ut eae us ly rn us ) (Ze 12a) 
pee ...(0) (0) (0) 

V = vy, Xx + Vo Y + V3 CZ .0-25) 

eo) k p() x x p§y A ps”? ee) 
a (0) (0) (0) 

oo). = ql x n a5” (2.12e) 


Substitution of these values into the zeroth-order equations 


and boundary conditions (2.9) gives the following solution 


Se 


y 9) = uw) = E/A (2.15a) 
gg QE a eg (2.13b) 
pool 3 app (2.13c) 
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i) = ES een (202d) 


(0) (0) 
Zs 


X/A +q =e Ae A COS. (2.13e) 


10 
ll 


(0) 


where the coefficient q5 was found from the additional con- 
dition that the bow shock is attached to the leading edge 


of the body. Recalling that the equations of the surface 


and the shock are given by S(x,y,t) = -yte(x-hcos?) = 0 
and G(x,y,t) = -yt+xtangd+eQ(x) = 0 we get the above result 
by letting x = 0 in the equation ¢(x-hcos8) = xeano+e (Q*-” 
mayo '*?). 
Similarly to solve the second boundary value 
problem we assume a solution of the form 
OEE es) (1) (1) 
U = uy Ux + us y + uz (2.14a) 
ce) (1) (1) (1) 
V ov, x t voy + Vv, (2.14b) 
a (1) (1) ele) 
P = Py xX +P, ‘y + P3 (2a) 
SE) es (1) (1) (ly) 
R poem er, OY + ©, (2.14d) 
g(t) = ann x + ast) x + ayy (2.14e) 


and substitute these values to the first-order equations and 


boundary conditions (2.10). The resulting solution is given 
below 
ee pa (Ll) ate i 
uy = [Cq, /M + (E +D/M,)/A + 2tand] 
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. = ta\7) (BHC/M ) + (E+F+D/M_ WAY cane 
ao = MeoS EEA =F) 
- = dl, cae = - hcos 6 
am = (aq!?) + Bya - 1)/tan 6 
im = ora + D/A + 2M, tan 9 
»” ~ 
_ =n cos 6[C(B-1)/A -D] 
a = i + (D+C-G)/A + 2 M. tan $ 
a = Noe -c)q\*) + (¢ +H -C -D)/A]/tan o- 2M. 
a = h cos 8 [1G (B-1) /A - H] 
a = be 2 (A-B-C+E/M_ \ JA -M. tan (D/A+2 M tan?) 1/2(k7a 
+ OM tan ¢) 
= = Mo/(Mo-1) 


2. Two-Dimensional Oscillating Flat Plate--Expansion 


Side 


This problem is formulated and solved in [Ref. 8]. 
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a. Problem Formulation 

Consider a two-dimensional flat plate of length 
2 in a supersonic/hypersonic uniform steady flow of an 
inviscid perfect gas with constant specific heats. Assume 
that the plate is performing a low amplitude and frequency 
harmonic oscillation about its leading edge. Let a system 
of Cartesian coordinates Oxy be attached to the body so that 
O coincides with the leading edge of the plate and axis. Ox 
is along itS mean position (Figure 2.2a). The bow shock is 
assumed to be attached to the leading edge and the flow 
quantities on the upper surface of the plate (in Region C) 
are required. 

me Method of Solution 

The unsteady flow, over the upper surface of 
the oscillating plate, will be found by perturbing the 
steady Prandtl-Meyer flow, considered in Section II.A.4.b. 
Denote by U,, p, and ep, the velocity, pressure and density 
in Region A. Denote by Uy, Py Py the velocity, pressure 
and density of the reference steady flow in region C. The 
solution procedure that follows is similar to the one pre- 
sented in Section III.A.1.b and most of the assumptions and 
results given there, apply to this section too, provided 
erat u , Pp 


fe) 
ties over the compression side of the flat plate) are replaced 


, o., M_ and a_ (the reference steady flow quanti- 
O O O O 


py’ Py? Py My and aby (the reference steady flow quantities 


Over theexpansion side of the flat plate). 


by u 
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Figure 2.2a. Oscillating flat plate--Expansion side 
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Figure 2.2b. Oscillating flat plate--Polar coordinates 
for Prandtl-Meyer flow 





Non-dimensional time and lengths are introduced, 
defined by (2.4). As a result of the oscillation of the 
plate assume that the perturbed flow quantities in region C 
are given by (2.5). The resulting perturbation equations are 
given by (2.6). Assume that the perturbation quantities are 
given by (2.7) and let the time independent quantities U, V, 
P and R be expressed as power seriesin (1k) by (2.8). The 
zeroth and first-order equations, derived as before, are 


given by (2.9a-2.9d) and (2.10a-2.10d) restated below. 


(0) _ _ 1 (0) 

J = -— PL QealGa) 

1 

Oo (0) 

vi . (2.16b) 
x x 

yl) 4 yl) gw pl?) 2 9 @uilea 
x y ae 
(1) nl) Pe es (0)) 

ut) + mm a (2.17a) 

yt) + _ p (1) = y $9) (217) 
Xx My 4 

pil) - pM) 2 pl _ 20) (2.17¢) 
x x 

uit) gy) gy pO) 2 Ly RM) (ney 
x 4 lea ak ik 


We now consider the boundary conditions applica- 


ble to the problem. 


qe 





Let the equation of the oscillating surface of 


the plate be given by S(x,y,t) = y-ex = 0.2 The flow 


tangency condition (2.2), with V = u,{l +eU,ev}, gives after 


linearization 
p AE y = 0 (mies) 
yt) = xX At y = 0 (2.17e) 


We next consider the boundary conditions across 
the surface that 1s separating regions B anc C. We assume 
that as a result of the small amplitude slow oscillations 
of the body the separating surface is slightly deformed and 
1ts equation is given by G = -y +x tano +eQ(x) = 0 where 
Q(x) is an unknown function which may be expanded as 


Cer eyo (2) 


Q(x) =Q + ... . We call the flow expansion 
an expansion front or, simply, front and assume 

that upstream of it the Prandtl-Meyer flow is not disturbed 
while, along the front, the unsteady flow matches the steady 
Prandtl Meyer flow continuously. The assumptions made are 
completely analogous to the assumptions made in the case of 
a finite compression shock discontinuity and the Rankine- 


Hugoniot conditions (2.3) may be used to give the boundary 


conditions across the expansion front. We note that, since 


ithe difference in the form of this equation and the one 
considered in the previous section is due to our assumption 
that in this case the flat plate is oscillating about its 
LE and thus h = 0. 


Ue 





we are dealing with an infinitely weak discontinuity, the 


flow across the front is isentropic. The procedure is given in 


Appendix C and the resulting boundary conditions, after lineari- 


zation and use of equations (2.8), are given below. 


v9 = arg: (9) At y = x tand (2.16£) 
p69) 2 erg: (0) (2.16g) 
Te ee eg Oo (2.16h) 
Rie = Gigt NSE) 
ha = A'Q' (1) + B1Q 9} At y=xtan 6 G2) 
pit) = crgt tl 4 pg) (2.179) 
te = evo) + Fig?) (2.17h) 
eeemeetou =) + 3799) (2.174) 


where the coefficients A' through J' depend on the reference 
steady flow and are given in Appendix Cc. 
c. Complete Solution 
Two boundary value problems have, again, been 
set up and will be solved successively. The zeroth-order 
equations and boundary conditions (2.16) constitute the 


first problem, which corresponds to a steady 
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Prandtl-Meyer flow problem. To solve it assume, as before, 
that the unknowns have the form given by (2.12), substitute 


in (2.16) and get with K = M,/ Vue -1 


(0) (0) 


P = Pz = kK (2.18a) 
R69) = ane = fe (Cn lL siey), 
aoa) Ss Sem (2.18¢) 
3 Ll 
oi ree ial (2.184) 
gi? = is Sigs eet (2.18e) 
gq 
where the coefficient i =O Sance- ain this case, the 


plate is oscillating about its LE. 

The second boundary value problem consists of 
the first-order equations and boundary conditions given by 
(2.17). In this case the unknowns are expressed by (2.14) 


and the resulting solution is 


2 
M. (M2-2) 
pie) = St x - amy (2.19a) 


( Vm?-1) 3 


IL 


i nee 5 (1) (2.19b) 


2 
Vb) L _— — y (2519¢) 
1 


i dyOaaye M 
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y (27.9G)) 


3. Three-Dimensional Oscillating Wings of Arbitrary 
Planform Shape 


The problem is formulated and solved in [Ref. 5] 
and only the major steps in the method of solution are given 
ifere. 

a. Problem Formulation 

Consider an oscillating wing of arbitrary plan- 

form shape at an angle of attack ina steady, uniform, super- 
sonic/hypersonic flow (Figure 2.3). Assume that the oscilla- 
tions are periodic with small amplitude and frequency and 
that the bow shock is attached to the wing. We let the 
pressure, density, velocity and Mach number of the approach- 
ing flow be given by p., 9,, UL, and M_. We also denote the 
total area and the root chord of the wing by S and % respec- 
tively and the distance of the pivot position from the leading 
edge by Xo: We assume that the pitching motion of the wing 


is described by 


Q(t) = @6 sue 


where 6 and w are the amplitude and circular frequency of 
oscillation and t is the non-dimensional time. 
We define the reduced frequency of oscillation 


by 


k = wk /U 


UE 





ty 





x/I| = f(y/b) ————— 
i x/l=g(y/b) 
Uco x 
oe 
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Figure 2.3. Three-dimensional wing 
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Dee Metnod of Solution 
The in-pitch stability derivatives of the wing 
are required. The pitching moment coefficient Ca! the 
stiffness derivative =e and the damping in-pitch deriva- 


Q 
tive aC... are defined by 


2 
eee = 2M/pUL2S (1/28) J) (ax) 6, (try t) as 


1 
<D 
ct 

d 

) 
+ 
}- 
a 

! 

‘o) 


CZ 732.0) 


where M is the moment about the pivot axis and oe is the 
pressure coefficient. 


The pressure coefficient is defined as usual by 


2(p -p,) 
Cc = 
aes 


and using the two-dimensional flow assumption can also be written in the form 


: x 
co . 2S Cc 
a = ee ID Cc +)! 
= Pp’. + eee 7 (Bx Cx.) ] G22) 


where A, B and C are dimensionless constants which are 
functions of geometry and steady flow quantities and (C6 
is the mean pressure coefficient corresponding to @(t) = 0. 
This coefficient has no contribution to the stability 
derivatives and will be neglected. 


Our goal is to calculate the coefficients A, B, 


C appearing in relation (2.21) and use the pressure coefficient 


we 





m@merelation (2.20) to obtain the stability derivatives. 
It is clear that, for a two-dimensional flat plate with 
attached shock wave, the coefficients will be the sum of 
independently obtained coefficients over the compression 
and expansion sides of the plate respectively. Thus the 


coefficients will have the form 


A = A, -~ A, Pena jmce= sc = CU, (Ze 2) 


where the subscripts 2% and u denote the lower and upper 
surface respectively. 
The results from Section III.A.1 and [Ref. 6,7] 


are used to obtain coefficients Ay, By and Co of the form 


~~ te 


Ay = AC7Ar By = Wo (2G-TI), Cy = uot C2 2 ok} 


where a and He are coefficients introduced to account for 


differences in notation given by 


_ ea ay 7 
r = M5, Zz , Me = d Uan/ 4, 


and G, I are quantities defined in [Ref. 6,7]. 
Similarly the results from Section III.A.2 and 


[Ref. 8] are used to obtain coefficients AL Be and Se Git 


the form 
Pee = Me 7-1) °°? Cee (2.23b) 
u ge 1 j u oh Le Ib 
7 2 = 
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where hye H, are correction coefficients to match present 
notation of the form 
- = om Gh? , Hy = AYU /u, 
ae scOlution 
BuoSeeemcang Lelations (2.23) in relation (2.22) 

and putting the resulting expressions in relation (2.21) 
we Obtain the pressure coefficient for a two-dimensional 
flat plate. We shall employ the strip theory to solve the 
three-dimensional wing problem in hand, which means we will 
assume that the flow, at each point of the wing, is two- 
dimensional locally. This assumption permits the use of the 
two-dimensional flat plate pressure coefficient for the case 
of the three-dimensional flat wing. 


Using relations (2.21) and (2.20) we get 


“Ch = ACT, - x,/2) (2.24a) 
2 
“C= (BI, +(C-B) Ig] -[ (BHC), +(C = B)I 51x /2+ C(x,/2) 
: (2.24b) 
where 
l l 
Mes & f (g°-£-)dn, 1, = 2k/3 f (g?-£°)an (2.25) 
l ; z 5 
1 ee 
T, = 2k, f{ f(q-f)dn, 1, = k Jf £(g°-£°)dn 
3 aa 4 5 
k = b/s, fae = Yb 
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For delta wings with power law leading edges, i.e., wings 


with f(n) = nt7n and g(n) = 1, equations (2.25) become 
I, = (nt+l)/(n+2), I, = (n+1)/(n+3) (25254) 
I, = n/(n+2) , I, = n/(n+3) 


The minimum value of Ca, Obtained for pivot position 
8 


x 7k = [(B+(2n+1)C]}/2C (n+2) 


is 


(-c =~ (Beem DC 2 4 4 (Btnc)/(n+3) 
j 


By setting (-C_ )_.. = 0a stability boundary for power law 
ma min 

delta wings may be obtained, which is practically independent 

of the power n, as shown in (Figure 2.4). Plots of Ch, and 


9 
G vs a for several values of pivot axis position and powers 


m 
n .. eawen in (Figure 2.5). In (Figures 2.6, 2.7) the 
Stability derivatives vs the pivot position are plotted for 
several values of power n and angles of attack 10° and 20°. 
Comparisons of results obtained by the present theory with 
results obtained by other theories and related discussions 
are included in [Ref. 5] and will not be repeated here. 


Nevertheless in Section III.C we will compare the present 


results with potential flow theory results. 
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Figure 2.4. Stability boundary 
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Figure 2.5. Stability derivatives vs AOA 
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Figure 2.6. Stiffness derivative vs pivot position 
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Be. LINEARIZED THEORY 

It was pointed out in Section II that several methods 
are used in the linearized theory of supersonic aerodynamics. 
A simple application of one of these methods, the so-called 
method of fundamental solutions, will be described in this 
subsection. For this application a distribution of pulsating 
sources over a "Simple" planform will be considered. 

Before dealing with this application, however, the notion 
of a disturbance propagation in supersonic flow will be re- 
viewed and the fundamental solution for a moving pulsating 
source will be introduced. 

ie ecopagation of Disturbances 

Supersonic flow is dominated by the fact that dis- 
turbances travel with finite velocities, namely, the speed of 
sound. In formulating the linearized potential equation 
(1.28) deviations of the speed of sound from its free stream 
value a, were neglected and also the perturbation velocities 
u', v', w' were taken very small compared to the free stream 
velocity U.. As a result any disturbance from a source 
located at point (x,y,z), in a coordinate system fixed to 
a body in the flow, can be felt only inside or on the surface 
of a right cirular cone whose axis points downstream from 
the source (Figure 2.8). An observer moving with the fluid 
sees a pulse emitted at t = 0, expanding on a spherical sur- 
face with instantaneous radius a t and center moving down- 


stream with velocity U,. The positions of these expanding, 
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Upstream zone of influence Downstream zme of influence 





Figure 2.8. Upstream and downstream zones of influence. 
Disturbance located at (&,n,t) at time t = 0 


moving spheres form an envelope which is a cone of semivertex 


angle 


u = sin ——- = Sin 


where yw is the Mach angle. 
This cone is known as the Mach cone or downstream 
zone of influence of the point (x,y,z) and its equation is 


given by 
(E=x)* - (M2-1) [(n-y)* + (¢-z)*] = 0 Os 


On the other hand the point (x,y,z) can be influenced 
by sources whose locus is evidently a similar cone directed 
forward from (x,y,z). This cone is known as the forecone 
or upstream zone of influence and has the same equation except 


Meat cor this case x > &. 


86 





Assume now that a steady supersonic flow has been 
established past a stationary three dimensional wing lying 
very close to the xy-plane. Regarding each point of the wing 
as a disturbance source we can see that the downstream zone 
of influence for the entire wing is bounded by the envelope 
of the Mach cones emanating from the leading edge. If the 
leading edge is straight (or the wing is two-dimensional) 
the envelope reduces to the so-called Mach wedges. 

To calculate the fluid motion at any point (x,y,2Z) 
we need to consider only the contribution from the disturbance 
sources that belong to the region of the xy-plane intercepted 
by the forecone from point (x,y,z). This area of influence 
forms a hyperbola and, with sources assumed to lie on the 
xy-plane, is found from the equation of the forecone (2.26) 


by setting 5 = 0. We thus get 


= yt V (x+)¥MS-1) -27 (2.27) 


ie 2 


A generalized supersonic planform with leading edge 
AA'C'C, trailing edge DD'F'F and streamwise tips AF and CD 
is considered next (Figure 2.9). For each point (x,y,0°) the 
area of the sources that influence the point reduces to a 
region bounded by two straight lines upstream of the point. 
These lines found from the forecone equation (2.26) by setting 


z= 0 and ¢ = 0 are given by 


E-x = ¢ VM*-1 ena SoS Ox 
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Figure 2.9. Generalized supersonic planform 


and are shown as dotted lines making an angle u with the 


maxis (Figure 2.9). 


The portions A'BC' and D'EF' of the leading and 


trailing edges are called supersonic since the velocity normal 
to these edges is greater than the speed of sound. Similarly 
the remaining portions of the leading and trailing edges are 
called subsonic since the normal component of velocity 1s 

less than the speed of sound. Along the supersonic portions 
of the leading and trailing edges there is no communication 


between the upper and lower surfaces of the wing and the flow 
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over the top or bottom of the planform can be calculated 
without reference to the shape of the opposite side. On 
the contrary the upper and lower surfaces of the wing are 
not independent along subsonic portions of the edges. 

Points 1, 2 and 3 Shown (Figure 2.9) are selected 
to illustrate different sorts of upstream influence regions. 
The difficulty in solving the linearized flow problem, i.e., 
finding the fluid motion at these points, increases as we 
meve £rom point 1 to point 3. 

There exists no universal method of approach in solving 
the linearized problem for different sorts of influence re- 
gions. Thus each planform shape calls for a different method 
Of approach. 

It is the simplest case of a planform with purely 
Supersonic leading and trailing edges, the so-called simple 
planform, that will be considered later in this section. 

2. Fundamental Solution of a Moving Source 
For U_ = 0 the linearized potential equation (1.28) 


reduces to 


CE 


y + Y + Y = - ty 
XX yY Cag a 


This is the classical wave equation for the propagation of 


sound in a still medium and its fundamental solution is given 


by 


¥(x,y,Z,t) = = F(t -=) (2.28) 
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where a. is the speed of sound of the undisturbed flow and 

me 1S the radial distance from the origin, i.e., r = Vx-+y* +2-. 
A solution of the linearized potential equation can 

thus be readily found if we can transform it into the classi- 

cal wave equation. To achieve this the following transfor- 


mation, known as the Lorentz transformation, is employed. 


~*~, 2 =1y, Gao Me 


Ss ae x) 
a (1 -M%) 
The velocity potential of a sound source fixed with 


respect to the xyz system is then found to be 


ua V1 =m elle xe 


nae, 7 ,Z,t) = ater ee f(t as ——_—_—___—_——_—_, 


Lee = 2 


el Mi.) 
Jt =u 
where R = Vx? + (1-M4) (y* + 2°) and 1-m? is known as the 


meanatl factor. 

To find the constants a and uw introduced by the 
Lorentz transformation we require that the sound source should 
produce constant sound flux independent of the free stream 


Mach number and we get 


a = a (l-M2) , u = Vi-m? 


Thus the velocity potential of the moving source becomes 


if in addition the relation A= - 7 is used, 
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iL piled aa ds 
rx ,;Y,2,c) = = aR er er gar 


Expressing the above result in a coordinate system moving 


uniformly with the source located at (&,n,2t) we obtain 


1 Bn an 1 
Tre Pit - > = id Ta F(t-T,) C2%.29) 


co 


ey y,z, t) = 


where: 
R = V(x-6)*+(1-M*) E(y -n)*? + (2-2) 71 


and 


oe a) te 
° a (1 M2) 
By comparison of expressions (2.28) and (2.29) it is seen 
that the solution for a moving source can be obtained from 
the solution for a stationary source by replacing the ordinary 
distance r by R in the amplitude and by D in the phase. The 
quantities R and D are called amplitude and phase radii 
respectively. 

A geometric interpretationis given in Figure 2.10. 
At time t a field point Q and a source O moving with velocity 
om in the negative x-direction are considered. For supersonic 
flow there are two spherical waves passing through point Q 
at time t. These waves originated from the source at times 


t-T and t-T at which times the source was located at 


QF 


and P. shown (Figure 2.10a). For subsonic flow 


positions P, 5 


oe: 





Keo 7) ee 

oe 
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Figure 2.10a. Supersonic source 





Pepiees 2-105. Subsonic source 


a2 





there is one spherical wave passing through point Q that 
Originated at time t-T when the source was located at point 


P (Figure 2.10b). From the geometry of the figure we get 
(x-E-UT)* + (y-n)* + (2-5)% = aft (2.30) 


Solving this equation for time T we get for subsonic flow 
One real positive solution and for supersonic flow two real 
positive solutions. 

Physically we are looking for the effect that a dis- 
turbance, originating at point (&,n,f) at some time t-T, 
will have at some later time t at a point (x,y,z). In this 
sense the potential is a retarded potential. For supersonic 
flow the disturbance is first felt at some point (x,y,Z) 


after a certain time T, has elapsed. The point (x,y,2Z) 


| 
penetrates the wave front of the disturbed region and be- 
cause, it 1S moving at a speed greater than that of the wave 
front it emerges from the disturbed region at some later time 


iT For subsonic flow once the point (x,y,z) penetrates the 


2° 
wave front (after a certain time T has elapsed) it will remain 
in the disturbed region since its speed is less than that of 
the wave front. Finally the nonexistence of positive roots 

Ot Pe uation (2.30) should be associated with an undisturbed 
region, i.e., with ¥ = 0. 


In view of above reasoning the source solution for 


supersonic flow takes the form 


9.3 





¥(x,y,z,t) = - aaRlF (t-T,) +F(t-7,)1 


Wwnere: 
D 1 Ma (*-6) -R 
aoe oo. 
00 00 M“-1 
D 1 Mae) ok 
mec a. a! 
o Me-1 


R = V(x-€)* = (Mo-1) [(y-n)* + (2-2) 71 


For a purely harmonic time dependence the source 


solution takes the following final form 


oS ame 5 
mm ,Y,2,¢) = = san e cos 5 R (2534) 


where: 


W is the frequency of oscillation 


wM 
y= ——— is the compressible reduced frequency of 
a (Mo-1) oscillation 
R = V(x-e)* = (Mo-1) [(y-n) * + (2-2) 7] 


a Simple Planform Solution 


Consider a three dimensional wing performing a small 
amplitude harmonic oscillation of circular frequency w. Let 
the wing surface be very close to the xy-plane of a Cartesian 
coordinate system attached to it and let the body move with 


Supersonic velocity U, in the negative x-direction. 
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The equations of the upper and lower wing surfaces 


are given by 


2 Ot 


S|, (*-y,t) Z-Z (x,y,t) z-h (x,ye 


Wt 


Th 
S, (x,y,t) Z ~Z, (x,y,t) z-h, (x,y)e 


and the linearized tangency conditions given by (1.29) become 


iy .ot.t) = eft tion +0. (2.32a) 

aaa me TBE a ae 
| 3h 

w(x,y,0,t) = e 9" [inh +U, =I (2.32b) 


Assume for simplicity that we are dealing with a 
Simple planform and regard each point of the wing as a pul- 
Sating source. Recalling definitions and terms used in 
Section III.B.1 this means that 

a) The leading and trailing edges of the planform are 
purely supersonic and the flows over the upper and lower 
sides of the wing are independent. 

b) Finding the fluid motion at each point of the surfaces 
involves the same sort of upstream influence zone, namely, the 
Sere indicated by point l in Figure 2.9. 

c) For a general point (x,y,z) the contribution of the 
disturbance sources that lie on the hyperbola with end points 
ny and No given by equation (2.27) is to be considered. 

In view of the above the complete solution for the 


velocity potential can be found by superimposing the effects 
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of all sources that lie within the upstream influence Zone 
(Figure 2.11). Integrating the fundamental solution for a 


moving source, given by relation (2.31), over the hyperbola 





we get 
— <> oe ) 2 — 
aie BRE Vd nD | nye rw (x-§) 
¥(x,y,2,t) = ore f f Be 
E=0 ny 
Cos M_ R 
RB qadé& dn 225.5) 


where q(&,n) is the unknown source strength. 





simple planform 
f orecone 





(a) (b) 


Eagurer2.l0. Limits of integration 
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To find the source strength we note that in the pro- 
cedure followed the wing was represented by a pulsating 
panel embedded in an otherwise rigid xy-plane. A mathemati- 
cal boundary value problem has thus been formulated for which 
the normal velocity of the fluid on the panel is prescribed 
by relations (2.32) and outside the panel it is zero, since 
no pulsating sources exist there. 

It can be shown based on physical considerations that 
the source strength can be expressed in terms of the normal 


velocity by the realtion 


a(gsn) = 2w(E,n) = 2(iwh +U, 2) 


Substitution of this value in equation (2.33) per- 
mits a straightforward calculation of the velocity potential 


Ee 


C. COMPARISONS WITH LINEARIZED POTENTIAL FLOW THEORY 

In this section we consider sweptback tapered wings with 
straight supersonic leading and trailing edges and streamwise 
tips (Figure 2.12). For this planform shape equations (2.25) 


become 


; De 

I, = 4B (k, ki)/3 +k, +11/8 (2.34a) 

I. = 228[(ko-k>) /4 +k*> +3k./2 +11/38 (2.34b) 
2 a wll 2 2 

I, = 228k, [(k,-k,)/3+0.51/s (2.34¢) 


Sil 





f(n)=k,n 


Cy 


(x5, 6) 


(1 ,O) 


q(n)=k.n vel 


X 


Figure 2.12. Geometry of sweptback tapered wings 


a Zee 
a &Bk, IL (kK, kj) /4 + 2k5/3 +0.51/S (2.34d) 
S/2£8 = k,-k +2, Ky = K/h, kK, = x5/k-1 

Some special cases readily follow. 

For a two-dimensional flat plate 
I, = AEF I, = Sy, it hegre OF S/2XB = 2 


and the stability derivatives become 


ee 
9 


A(0.5 -x /%) 


ec 
9 


Hy 
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B/3 - (BHC) (x,/2)/2 + C(x /2) é 





For the case of a triangular delta wing, 


iT eae / 3, iL =e. 1 ai By at Ee Ae o/i8 =. 1 


and the stability derivatives are 


j 
Q 
iI 


A(2/3 -x/%) 


i 
2) 
i 


(B+C)/4 - (B+3C)(x_/2)/3 + C(x /2) 4 


For oscillations at small angle of attack the coeffi- 


cients A, B, C appearing in equation (2.21) become 


~ = C = a B = a (87-1) /B (22.55) 





where 8 = /m?-1 ‘ana OS Wey eve eyocGineneieel Halim, Cenohantd= 
slope (00 = 4/8). 

In this limiting case the stability derivatives given 
above for a two-dimensional flat plate and a triangular delta 
wing become, as expected, identical to well known formulas 
based on potential flow theory (Ref. 10: pp. 52,144]. Asa 
result of the strip theory approximation the stability deriva- 
tives of a rectangular wing are independent of its aspect 
ratio and the stability boundary is always the one for a 
two-dimensional wing [Ref. 10; Fig. (RTA 

Formulas (2.34) and (2.35), based on Hui's theory, 
were used to calculate the stability derivatives for several 


sweptback wings, with straight leading and trailing edges 


oe) 





meena the results were compared with results based on potential 
flow theory, as given in [Ref. 11,12]. The comparisons are 
shown in (Figures 2.13-2.17) and the procedure followed to 
read values for the stability derivatives from [Ref. ll, 12] 
is described in Appendix E. Good agreement is generally 


shown. 
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Figure 2.13. Stability derivatives vs Mach number 
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Figure 2.14. Stability derivatives vs leading edge 
sweep angle 


HEORZ 








Ht oS. | THEORY 


TW AGCAN TN 2699. 3196 





0 f25 .5O 75 1.0 


Hines 2.Lo. Stability derivatives vs taper ratio 


iOS 





dae eo HUIS «THEORY 
_. + — + -NACA TN 2699, 3196 


6 Mz2, 7-25, L, p= 455 x/1 = 0 





Figure 2.16. Stability derivatives vs aspect ratio 
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Figure 2.17a. Stiffness derivative vs Mach number 
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IV. EXTENSION OF PROFESSOR HUI'S THEORY-- 
UNSTEADY , NONUNIFORM UPSTREAM FLOW 

Consider a uniform, spacewise and timewise, supersonic/ 
hypersonic flow of an inviscid, perfect gas, with constant 
specific heats, past a formation of two wedges at design 
Sondition (Figure SoA) oe The wedges are assumed to be 
oscillating in pitch, with small amplitude and frequency, 
in an independent fashion and the bow shocks are assumed to 
be attached to both bodies at all times. The flow quantities 
over the upper and lower sides of the second two-dimensional 
wedge (in Regions C and D) are required. It is clear that, 
depending on the difference of the semi-vertex angles of the 
two wedges, either 

a) both sides of the second wedge are compression sides 
weegure 3.la) or, 
b) one side of the second wedge is a compression side and 

the other side is an expansion side (Figure Soi o 

These compression and expansion side problems are con- 
Sidered in subsequent subsections. 

In subsections A and B the flow quantities over the com- 
pression and expansion sides of an oscillating wedge will be 


given. The upstream flow (in Region B), is assumed to be 


this problem, suggested by Prof. Platzer, is of interest 
in high-speed turbomachinery aerodynamics. 
5 : 


meee oe. possibility of the second side being neither a 
compression nor an-expansion side is not excluded. In this 
case though region C becomes an extension of region 8,in 
which region the flow quantities are assumed to be known. 
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Fig. 3.1a FOrmation of two wedges. Both sides. of 


wedge Li are compression. sides. 
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‘oscillatory with circular frequency w and flow quantities 


of the form 
= uo (1 + €U) Cele 
2 V CS) 
p = Pott + eyMoP) (ri lc) 
on p61 + eM,R) (31d) 
where: 
u, v are the x- and y-velocity components relative 
to the coordinate system shown (Figure 3.la); 
Ee = aoe with ¢€: a small parameter characteristic of 


the deviation of the upstream flow quantities 
from their average value; 


k: the so-called reduced frequency defined by 
wk 


k=—; 
u 
O 


@: a characteristic length of the flow; 


Us ePo Po and Mo are the mean or average velocity, pressure, 
density and Mach number of the upstream flow; 


U,V,P and R: are time independent quantities considered to 
be known functions of the non-dimensional 
Spatial coordinates x and y. 


We will restrictthe analysis to small reduced frequencies 


k and will assume that these quantities have the form 


G 
Il 


ye cixyull? = a9) 4 (ik) (aft) xtul?? yeu$*?) (moan) 


(1) 


eV 
Ih 


ee ve cixyv')) = vi) 4 (ik) Cv Y) yey St?) (3. 2b) 
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0 ali 0 
i = + ) 4 (ik) P| ) = Ps + (ik) (pit? xp!) yep?) co 2C) 
me Rs ik)Ri)) q ry) + (ik) (rp?) xer$)) ytr (2) ) (3. 2a) 
where the coefficients oe rar are known constants. 


It should be noted that the above forms contain as 
Special cases the closed form solutions found in Sections 
TIIT.A.1 and III.A.2 for the flow over the upper and lower 
Sides of an oscillating flat plate at an arbitrary angle of 
attack. Thus the results given in the following subsections 
A and B hold true for upstream oscillatory fields of the 
general form (3.2) and are not limited to oscillatory fields 
created by oscillating flat plates/wedges. Similarly the 
solution given in subsection A contains as a special case 
the solution to the problem of a stationary wedge in an 
Oscillatory, uniform spacewise, hypersonic free stream, which 
was studied in [Ref. 9]. Furthermore the solution presented 
here is exact and holds for the complete supersonic/hypersonic 
speed range. 

In subsection C an alternate approach to the expansion 
Side problem is suggested. In this approach the upstream 
flow is assumed to be oscillatory and, for simplicity, uni- 
form spacewise. Two boundary value problems with linear 
equations and boundary conditions are formulated and a closed 


form solution of the problem is sought. 
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A. UNSTEADY, NON-UNIFORM FLOW PAST AN OSCILLATING TWO- 
DIMENSIONAL WEDGE--COMPRESSION SIDE 


Consider a wedge (Body I), at design condition, osci- 
llating with small amplitude and frequency in supersonic/ 
hypersonic flow (Figure 3.2). This flow problem was studied 
in Section III.A.1 and the flow field quantities Uns Vas Pp: 
Pps ... in region B (expressed relative to axes O'x'y') were 
completely determined. Assume now that a second wedge 
(Body II), located entirely in Region B, is oscillating also 
with small amplitude and frequency. The flow field quanti- 


Eies Ua, Vor Pur Py, ++. in Region C are required.+ 
To solve the problem we adopt the following procedure. 

a) Express all flow field quantities relative to the 
coordinate axes Oxy attached to body II in its mean position. 

b) Assume that both bodies I and II are stationary and 
find the reference steady flow quantities in Region C. 

c) Assume that body I is oscillating while body II is 
kept stationary and superimpose to the steady flow quanti- 
ties in region C perturbation quantities due to the oscilla- 
mon of body I. 

ad) Assume next that body II also is oscillating and 
Superimpose to the aleady perturbed flow quantities in Region 


C new perturbation quantities due to the oscillation of body 


iL 


1 Plow quantities in region D can be found ina completely 
analogous way (by simply letting 0, = 8, +8, in the solution) 
Since the flows are assumed independent. 


IGE 





_ a . J a 
3 ot — 
o 
: : : . ; 
ae ee eee ee eee 
° 


ae of 


Nari ge5.2 


I Llustr ation 


HES: 


o f 


the 





flow quantities. 





The procedure described above effectively breaks down 
the flow problem in hand (namely, the calculation of the 
flow quantities in Region C when both wedges are oscillating) 
into a sequence of three separate problems which may be solved 
successively to give the final results. 

The solution method, which is based on Professor Hui's 
theory, is presented in Sections IV.A.2 and IV.A.3. The 
final results are given in Section IV.A.4. In the following 
section IV.A.1, the problem is formulated and the flow field 
quantities in Region B are expressed relative to axes Oxy 


attached to body II in its mean position. 


imme ormulation of Problem--Expressing Upstream Quantities 
Relative to Coordinate System Attached to Second Body 


Consider steady uniform supersonic/hypersonic flow 
past the formation of the wedges described above (Figure 3.3). 
Assume that the wedges have chord lengths 24 and Lo and are 
oscillating with circular frequencies Wy and Wo Cartesian 
coordinate systems O'x'y' and Oxy are attached to the bodies 
with origins placed at the wedge apexes and axes O'x', Ox 
along the mean positions of the upper surfaces. Let the 
steady flow quantities in Regions B and C (stationary wedges) 


be given by u Po 4 Py and sie Os. Por respectively. Also 


’ 
ore “3 B O Oo 


denote by t the time variable and by u,, Ver Pp and Pa the 
perturbed quantities in Region B expressed relative to axes 
O'x'y'. The departure of the perturbed flow from the steady 
flow in Region B is characterized by the small quantity 55, 


introduced below and the reduced frequency associated with 
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the Beertlation of the first wedge is defined, as usual, by 





Te a 
Ky = ——. Similarly the reduced frequency associated with 
Oo TK 
the oscitlation of the second wedge is defined by kK, = < a 
O 


-Non-dimensional time and lengths are introduced associated 


with bodies I and II as follows, 
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We assume that the perturbed quantities in Region B 


have the form 


un = u (1 + €,U) (35a) 
B 
Vn = Ey Boe V (3285p) 
D. = Fa.‘ + ce (Sec) 
Pp. = op. (1 + €,M_ R) (3.3d) 
B On iL Op 
mn ik,t 
where Ey = e,¢ and 
ek) (ul?) x? eu ys tu?) (3.4a) 
3 ie iL 2 S 
a (0) : al te Ce ane (1) 
y= V3 + (ik,) (vy eT AP ie Pare ) (3. 4b) 
- (0) it) ae (ole) aes, iL) 
P = P, + (ik, ) (p, x +P, y +P, ) bake) 
iL 
R = + (ik,) (rit) x: try) y' rs - (3.4d) 
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with coefficients u --. assumed to be given. We 


v 


ae 
recall that in Sections III.A.1 and III.A.2 the problems 
of uniform, steady supersonic/hypersonic flow, past harmonically 
oscillating wedges/flat plates were studied and the perturbed 
quantities behind the shocks/expansion fans were, in both 
cases, found to be of the form assumed above by relations (3.3) 
and (3.4). Thus the analysis that follows does not distinguish 
whether the disturbed flow field in Region B has originated 
from an oscillating wedge or flat plate over its expansion 
and compression sides and the appropriate coefficients and 
parameters should be used to make the distinction for the 
specific case considered. 

Next we express the quantities given by (3.3) and 
(3.4) relative to the coordinate system O xy attached to 
body II. Parallel transformation and rotation of axes gives 
for the general quantity 


oe 


(ee) 
ie oe 2 


(le) 
3) 


me (0) ; 
M = mM. + (1k, ) (m m von 


where symbols M and m stand for capital and lower case sym- 


mos U, V, ... and u, v, .-. respectively. 
Maen” & GGk.)c" x tWer y+e! ) G5) 
S i m4 ma mM. 


where 


OS = Oy cos o4 +B, sin oh Bo = By cos O41 — a, sind), 


Pay 





d Z IL 
Cn = (my*) cos 6, + my) sin 6.) Ro/k, 
“no = (ms*? cos 94 -m)*) sin G4) 4/2, 
Sn = my) + (amy) +8 mi1) 7/7, 


Combining relations (3.3) and (3.5) we express the perturbed 


Guantities Bs and Pp relative to axes 0 xy in the following 


form 
Meee Dp ilt+e_(C_ +(ik,\C_ x+C_ ytc_ ))} once 
p Op sane bs By Ps 
9, = ep. {lt+e,[(C_ +(ik,)(C_ x#C_ yto. )]} (3.6b) 
B On 1 ly IP ry ro i 
with 
(0) ep) 
C = y M p ’ C = M E 
ES Op 3 ry oO; S) 
eC = ¥yM_C' , C = MC! (juaaily,.2 ; 3) 
Ps On P, rs On r. 


To find the perturbed velocity components Ups ve we take 


their components along the new axes and add them to get 


& 
ll 


5 oe 8, +e, (U cos 68,+Vsin9,)] 


< 


= u {[{-siné 
fe) 


: +e, (-U sin6,+V cos 64) ] 


d 
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and finally using relations (3.5) we have 


u, = u_ icos 6,+¢e,[C_ +(ik,)(C. x+#C_ y+C_)]} (3.6c) 
B On d i Us a uy u, us 
Vv, = u_{-sind,+e,[C_ +(ik,)(C_ xtc. ytc )]} (3.64) 
B On d Jk ve ik Vi V5 V3 
where: 
(0) 0) ae (0) =. (0) 
— - = -_ 
eo. u, cost. V3 sino, cy. u, Sing +v, cost, 
— q q . — 
7 = Cy, cose atc’ sind. (ay = 1, 25:3) 
3 J J 
as -c'' ] ' a 
PY. Ci, Sing s+c) cose. (5 peep oS) 


Relations (3.6) give the complete set of flow quantities in 
Region B expressed relative to coordinate system Oxy. 
feeeecrcurbed Equations of Motion 

In this section we derive the perturbed equations of 
motion in Region C, when both wedges are oscillating (Figure 
ez). 1f10 £L£ind the form of the perturbed quantities in this 
region, relative to axes O xy, we proceed as follows. 

First we assume that both bodies are stationary and 
denote the (reference) steady flow quantities in Region C 
by Yor Por P57 Mo (Figure 3.4a). Obviously veo OVror the 
coordinate system chosen. 

Next we assume that body I is oscillating while 
body II is kept stationary. Due to the oscillation of body 


I we superimpose to the field quantities perturbations of 


the form (Figure 3.4b), 
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9. 3.4bD Flow quantities in Region C, Body | 


oscillating while Body Il is stationary. 





Me = ce tty, - 
e1¥, = Pe te = 
€1P) 7 get Ap ume, 
E104 7 ye tok, 


€,U5U, (oa) 
EUV, (S272) 
EPS YMOP) (3G) 
EP MOR] eo. 7d) 


Finally we assume that both bodies are oscillating. 


Due to the oscillation of body II we superimpose to the field 


quantities perturbations of the form 


EU, = E5€ Ze aeU = 
By, = ¢,¢ 7 uv, - 
QP, 7 S302 pM ,P 
E5Py a Se 2 21H, 


The flow field quantities in 


become (terms of order EjEn 1 FE, £5 


u= uo (1 + €,U, + €5U.) 
v= ule Vy + E5V5) 
p = poli + yMo(e,P, +e 


Z 


(Figure 3.4c) 


€ uu, ocd) 
EUV, (Sie os} 
E5PLYMP 5 (3.06) 

= E50 MOR, (320d) 


Req Oneeceromscne last «case 
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bodies oscillating, 
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0 = poli ~ Mo (eR, +€5R5)] (3.9d) 


where the quantities Us, U., Vie [te ace £uNncEIONS Of the 
non-dimensional variables x and y. For small Ky and kK, 
we may expand these quantities aS power series in (ik) of 


the general form 


k= a + (ik. a ee 1) (3.10) 
with X denoting U, V, P and R. 

To simplify we assume that the characteristic lengths 
Ry and L., are equal and we put (3.9) into the governing 
equations of motion (2.1). Using (3.10) and equating the 
‘terms of the same order of (ik) in each of the resulting 


equations we obtain 


m0) 9 D) Se p (0) 


l-lx * ©2%2x Wige Ba ose UM ae 
my, + jv (e,Py0 + oP oy) /My (3.11b) 
e, (Pio? = RIO’) + (PAY? - Ro) = 0 mie) 
en, ty tM QR Rio?) + 6, (U3) +V Voy tMQR3y) =0 (3.114) 

ey(k us Up +4, ut) se (ik, yu (U5)? +050’) 
= -u le, (ik,) Pi?) +e, (ik, yp 1m, ea) 
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(ik) (ay vl ea vy) p40 (iy u, (vl) ay) 


B i. Our 2x 
= -u, fey (ik, ) Py)’ +65 (dk,)P 52) 1M, (3.12b) 
ey (ky) fu, (Pp) -R Ry) au (ep) Rt) | 
+ 6, (ik,)u, (Py) -R3°) 4p lt) Roe) = 0 (3.12c¢) 
e, (ik) Cu, MQRp HY, (ute) + HVT) MRD) 
+e, (ik,)u (MOR, +0) ay ou +M RA2)) = 0 (3.12¢) 


Bem@ations (3.11 and (3.12) are the zeroth and first-order 
equations for our problem. For the case of a single oscillating 
body they reduce to equations (2.9, 2.10) given in Section 
ite n-l, by letting e, = Oi; Uy = Vy = Py = Ry = 0 (body I 
considered missing) and dropping subscripts B and2 in re- 
sulting expressions. 

To solve for the sixteen unknowns contained in 
Systems (3.11) and (3.12) we may proceed in two steps. We 
May first find the eight unknowns associated with the oscilla- 
Elen Of body eo To do this we will assume that body II is 
Hot Oscillating (€, = 0) and solve successively the above 
systems of equations subject to appropriately formulated 
boundary conditions for this case. Next, we will find the 


eight unknowns associated with the oscillation of body res 


Unknowns with subscript l. 
-unknowns with subscript 2. 
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To do this we will assume that body II is oscillating while 

body I is present but stationary (Ee) = 0) and solve successively 
the resulting systems of equations subject to boundary con- 
ditions appropriate for this case. The second problem has 

been treated in Section II.A and we will not repeat the solu- 
tion for this case. We note though that the problem should 

be appropriately reformulated to account for the different 
direction and magnitude of the approaching steady flow 

velocity as indicated by the ratio of steady velocities 


aL /u., appearing in equations (3.12). 
B 


3. Boundary Conditions 
For both wedges oscillating, let the equation of the 


surface of body II be given by (Figures 3.4b, 3.4c), 


S = “~y t+ e,(x - fh cos 8.) = 0 


ji 


and the equation of the shock attached to body II be given 


by (Figures 3,4), 
G = -y + x tango + €4Q, (x) + E5Q4X) = 0 


where Q, x) and Qa. (x) are unknown quantities associated with 
the oscillations of bodies I and II respectively and are to 
be found as part of the solution. 

We will formulate the boundary conditions for the 
case of stationary wedge II (€. = 0). 

The flow tangency condition (2.2) with V = uj{l +e, 1! 
€1V,3 gives 
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0 
vy = 0 JN i aa eae O (S73 a) 


i 
v1 = 0 At y = 0 (372235) 
To find the boundary conditions across the shock we 
use the Rankine-Hugoniot conditions. The procedure is given 
in Appendix D and the resulting zeroth and first-order condi- 
tions, after linearization and use of equation (3.10), have 


the following general form 


(0) + 10) 
aa 


a ad + Ky At y = xtan@o (3.14a) 
Se thy PK, KY ae +Ky a 
XX VY x fo) 


At y =x tan? (3.14b) 


where Y stands for U, V, P and R and the coefficients Koy 


Ky ; Ky , ee. are Known constants (functions of geometry 
x XX 


and steady flow quantities) given in Appendix D. 

gee oOlution 

Two boundary value problems have been set up and 

will be solved successively. The zeroth-order equations (3.11) 
with Ee, = QO and boundary conditions (3.13a) and (3.14a) con- 
stitute the first problem, which will give the steady flow 
quantities behind the shock, for the case of stationary 
wedges in uniform, steady supersonic/hypersonic flow. The 


first-order equations (3.12) with Ee, = OQ and boundary condi- 


tions (3.13b) and (3.14b) constitute the second problem. Since 


Ia 7! 





both problems are linear we assume that their solutions 


are linear combinations of the non-dimensional coordinates x 


ama y. 


For the first problem we assume a solution of the 


form 


| 
Us = 


(0) 
val 


mo) _ 
Ry - 


(0) 
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m) _ 
Qy aa 


Substitution 


and boundary 
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Dee ig 
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aye + a5" 


of these values in the zeroth order equations 


conditions gives the following solution. 


a3 a Nur qq 
0 

- RRx qq.” 
Be " Xx qq” 
ce = - wu x 


Bee 


(3.15a) 


(32E55)) 


(22 5S) 


(315d) 


(3.15e) 





For the second problem we assume a solution of the 


form 


y 2) 
iL 


ey _ 
Vi = 
cy _ 
Ry = 


io) 
al 


ce) 
Q, = 


1 3 
vit) x + vary ss a 
pte ey + 

(1) (1) (1) 


Pp) X + Po Y + P3 


(1) 2 (1) (2) 


> 


Substitution in the first-order equations and boundary con- 


ditions gives the solution 
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ult) asceen | 
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Dye ry 4 2 
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pit) x + pit) 
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WAS, 
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ae. (dL). 
26 _ Rex 72 7 SS 
(1 IL 
uy — - (Ps /M. + Bp Ye YQ) 
1 AL 
ws y= (-uy 4 Rudy + Wy) tan ¢ 
ay (1) 
3 = RuxF2 Z Ky 
It 1 
ae tee (as) 
r5 = ( ry + Red “ W,) /tan d 
’ - Keds) + Kp 
i 
49 ~ KY/K 
u 2 
fe) Ss) hoy oy Wo os 
O ‘e) 
2 
- Wy ee oe o “Sie ] 
tan 9 Mo tan o 
Wy = K. - Ky weve We Ky th 
XX VY ‘e) 


with Y standing for U, V, P and R. 
B. UNSTEADY NON-UNIFORM FLOW PAST AN OSCILLATING TWO- 
DIMENSIONAL WEDGE--EXPANSION SIDE 
Consider a wedge (Body I), at design condition, oscillating 
with small amplitude and frequency, in supersonic/hypersonic, 
steady, uniform flow (Figure 3.5). This flow problem was 


studied in Section III.A.1 and the flow field quantities in 
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two wedge 
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the flow quantities 





for the ~° 


- expansion case shown. 





Region B, expressed relative to axes O'x'y' were completely 
determined. Assume now that a second wedge (Body II), 
located entirely in Region B, is also oscillating with 
small amplitude and frequency. Let the semi-vertex angles 
of the wedges be Ya and ¥5 respectively and assume that their 
difference ¥ 5 sy is negative. In this case the upper side 
of the second wedge (Region D) becomes an expansion side and 
the exact solution given in the previous subsection does not 
apply. 
To solve the problem we adopt, as before, the following 

procedure. 

a) We assume that both bodies are stationary and find 
the reference flow quantitzes inside the expansion fan 
(Region C) and over the upper surface of the wedge (Region 
D). 

b) We assume that body I is oscillating, while body II 
is kept stationary and superimpose to the steady flow quan- 
tities in Regions C and D perturbation quantities due to the 
oscillation of body I. 

c) We finally assume that body II is also oscillating and 
Superimpose to the already perturbed flow quantities in 
Region C new perturbation quantities due to the oscillation 


of body a The problem is thus solved in three successive 


1 bisturbances Gue to the oscillation of body II propagate 
along the characteristics and are not expected to influence 
Region C except in the neighborhood of line OC. 
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steps. In what follows we will restrict the analysis to 
the case of a stationary wedge in oscillatory flow. The 
solution can be extended to the case of an oscillating wedge 
by following the procedure described in Section III.A.2. 

In subsections 1 and 2 the problem is formulated and the 
method of solution is presented. This method calculates the 
flow quantities in the expansion fan, in a sweeping fashion 
from Region B to Region D, along rays 9 = constant, with 
8 increasing in small steps from oe tomo (M)) (Figure 3.6) 
The procedure involves the repeated application of two basic 
steps and is described in subsection 3. The solution is 
approximate and becomes more exact in the limit as the number 


Of iterations increases. 


Pees ormulation of Problem--Expressing Upstream Quantities 
Relative to Coordinate System Attached to Second Body 


Consider steady, uniform, supersonic/hypersonic flow 
past the formation of the wedges described above (Figures, 
3.5, 3.6). Assume that the wedges have chord lengths dy and 
Le and are oscillating with circular frequencies w, and w,. 
Cartesian coordinate systems O'x'y’ and Oxy are attached to 
the bodies with origins placed at the wedge apexes and axes 
O'x4Ox along the mean positions of the upper surfaces. Let 
the steady flow quantities in Region B be given, relative 


pemaxes O'x'y', by u.,p_,p.,M ..- and in Region D, 
Sa Gay TR El 


Bemteivye to axes Oxy, by u, p, Pp, M, .-. . Let also the 
steady quantities along a ray 9 = constant in the Prandtl- 


Meyer expansion fan be given, relative to axes Oxy, by 
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body Ii. 





Ugr Voy Par Pos M, ... where 


u = we cos ([-8) + Vo sin (I-86) 


V = Vee Sin ([-8) a cos (T-8@) 


where the angle [ and the r- and 6@- velocity components 
_ and Vy were defined in Section II.A.4.b. 

We assume, as in Section IV.A.1, that the perturbed 
quantities in Region B have the form given by equations (3.3) 
and (3.4) and note as before that this form contains as 
special cases the form of the perturbed quantities behind 
oscillating flat plates/wedges. To express the quantities 
given by (3.3) and (3.4) relative to the coordinate system 
Oxy attached to body II we use relations (3.6) with 04 = Yor’)> 

2. Method of Solution 

Assume that the change in flow direction over the 

corner (Figure 3.7a) is obtained in n steps (Figure 3./7b). 


Assume that the step changes in flow direction are all equal, 


1, 


and also note that each of the line segments Dd Dy ay (oe=. ls, 
...-,n-l) extends to infinity. In the second case the single 
expansion fan has been replaced by n smaller expansion fans 


and the two problems are, physically, completely equivalent. 
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Next we assume that the change in flow direction over 
the step corners is achieved through expansion discontinui- 
ties, similar in nature to hypothetical expansion shocks, 
which we will call expansion fronts or, simply, fronts. 
Suitably spaced lines OB 


OB heey OB in (Figure 3.7a) 


is Ze 
correspond to the mean positions of the fronts Fi, Foy ee 
FO emanating from corners Di, Do, ae De respectively. 

These lines divide Region C in subregions Boe Bir Bor eee, 


Bo-1' By = D- The first subregion BA 1s separated from 


Nn 


Region B by line OB. = OB and from Region B, by front F 


al ig 


Subregions B, (k = 1,2,...,n-1) correspond to line segments 


k 
Dy) Dray and are separated from their adjacent subregions from 


the left by fronts Py and from the right by fronts F The 


k+1° 
last subregion BA = D corresponds to the line segment 
Del = DD (actual upper surface of wedge) and is separated 
from Region C by front Fae The flow quantities in this last 
Subregion are required. We intend to solve for the unknown 
quantities in this region by first obtaining the unknown 
quantities in Regions Boy Bor sees B-1 Successively. 

We observe that the original flow problem has been 
replaced by a set of n identical flow problems each of which 
involves the same step change in the direction of flow. We 
may further observe that the assumption of "jump" step changes 
in the direction of flow, via the so-called fronts, is sSimi- 
lar to the "jump" step changes in the direction of flow via 
ordinary compression shocks. Thus the two problems shown in 


(Figures 3.8a and 3.8b) are mathematically identical if the 
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Figure 3.8. Duality of problems 
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common term discontinuity is used to denote "jump" changes 
maeiLow direction and quantities across it. 

In Section IV.A the solution to the flow problem 
shown in (Figure 3.8a) was obtained. It was assumed there 
that the approaching flow was oscillatory of the form given 
by equations (3.3) and (3.4) and the flow quantities in Region 


BL were found relative to the coordinate system D We 


k*kYk* 
will see that the approaching flow for the dual problem shown 
in (Figure 3.8b) is of the same form and thus the solution 
given in Section IV.A.4 applies to the dual problem also, 
provided that the appropriate geometrical data and steady 
flow quantities are introduced for this case. The quantities 
found in Region By will be expressed relative to the next 
coordinate system De e414] and will be used as input 
quantities to an identical problem to give the solution for 
the flow quantities in the next Region Bead: The procedure 
can be repeated till the required quantities in the last 
region Dd, = D are found. 

Physically the assumption of a "jump" expansion dis- 
continuity is not accepted since such a finite discontinuity 
would lead, as explained in Section II.A.3.b to a decrease 
in entropy. The procedure described above, however, should 
give the exact solution in the limit, as the number of step 
changes in angle is increased. 

3. Solution Procedure 
Consider an arbitrary Region 3. (cele en) Deh ind 


front FL (Figures 3,8b, 3.9). Assume that both bodies are 


139 





B,=8B 
Q 
REGION 8B oS S (FAN ANGLE ) 
~ 
Re 
, * 


Pagure 3.9. 


Illustration of geometrical data 


140 





stationary and let the reference steady flow quantities in 
this region be Ur Ver Pur Pye Ms ~+. » These quantities 
can be found from the Prandtl-Meyer relations given in 

Section II.A.4.b. Expressed relative to coordinate system 


Ox, Y, these quantities are 


= Ve Oy) 

ae > Pz [E(M. ) /E(M,) | (3.17a) 
B B 

2 WA) 
a eee ge (Gers aie.) 
a = es, cos Ek + a Sin Sk (321 7e) 
i in = Sin SE ~ Yo. cos Ek (37d) 
M, = [(tan ne, /d)* #1] °°? (eee) 


where Ps ’ 25 ’ M. pa Mo Peneersligl= the reference steady FLOW 
B B B 
quantities in Region B and 


a pe’, E(M) = 1+ (y-1)M2/2 
Es 6 = _ = Il Deer O 5 
8, = a, +k at 8, = ae = tan ae. iL) JA 
= Pepe = 2 (n-k) a= ¥i-¥ 
ae kok! a n ’ a, 
es = c¢ sin MOL Vo = jr ¢ cos Gy 
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Assume that body I starts oscillating and denote the 


perturbed 


quantities in Region B. expressed relative to 


coordinate system Ox, Y,. by 


uw, = 


= 
rT 


where the 


R, may be 


oscillations, as power series in (ik,) of the form 


Yi = Y 


u, (1 - €,U,) (35,10) 
uy (€, Vp) (3.185) 
Pp, fl + €, YM, Py) (Sole, 
Pp, (1 a8 €, MR) (33a) 
unknown time-independent quantities U/, Vie Pye 
expressed, for small amplitude and frequency of 
1. 


, (0) 


. gatas) 
” ” + (ik,)¥y | eae 


ferme, Standing for U', V', P', R'. 


Assume that, with body I oscillating, the equation 


Sreene front FL 1s given by 


O(F,) 


where 


2 Ne tan >, +€,Q, (x,) = 0 (S29) 


k, is the reduced frequency parameter associated with 
the ostillations of body I and should not be confused with 
the integer k. 
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ee 2 Gt) > fo-o)k/n, = T = a+5+P(M ) 


[A (Me Sy - tan 


B B °p 
Assume that the approaching flow from Region Bry 
has the following form, with flow quantities expressed rela- 


tive to system OX, Yy1 
u_) = U,_p (cos 6, + €,U,_)) (3.20a) 
Yo = U,_z(-sin 6, + €,V,_)) @rZ0b) 
mes pay (i + €,R,_7) (3.20c) 
ome = «2, (1 + €,P,_,) (3. 20d) 
where O5 = (¥5-¥)/n and the time independent quantities 


Uriy7 Vyiyr Ruy and P,_, are assumed to be known quantities 


of the form 


Yeo] = Cy + (2k,) (C, x, + Cy Y, +Cy ) (322)) 
Oo 1 2 3 


Poseaey Standing for U, V, R and P. 

We observe that the problem formulated above for the 
arbitrary Region By 1s mathematically equivalent to the 
problem considered in Section IV.A since the approaching 
flow has the same form, the change in flow quantities occurs 
through a discontinuity that satisfies the same equation and 


therefore the same Rankine-Hugoniot boundary conditions, the 
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unknown perturbed flow quantities along the surface have the 
same form and should satisfy the same flow tangency boundary 
condition and the same governing equations of motion. We 
may therefore borrow from Section IV.A.4 the solution for our 
problem as given by equations (3.15) and (3.16). As a final 
step we should show that this solution expressed relative to 
the coordinate system OX. 1 YRal (associated with the problem 
that is to be considered next in order to find the flow quan- 
tities in the adjacent Region B41) is of the form assumed 

by relations (3.20) and (3.21). We see that this is indeed 
the case since the solution given by equations (3.22) below 
is expressed relative to the rotated system of axes OX, 1 YRel 
by equations (3.23). 


The solution given by equations (3.9), (3.15) and 


(3.16) is 
Be = u{l+e,(ur'”) + (At pag 9 5 (eee 2a) 
¥ = u,e,(vi'?) + (ik vil) (3.22b) 
Py = py {1 -- e,yM, (PY (°? + CCF: | Saleh E (3222¢) 
c= 0, {1 + em, (n°? + (ax, )me'*) 3} (3.224) 
with 
gy (9 _ ne yO) 2 9, pe) 2 (0), 
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moO) (0) BCE -0¢ (ele (Ir) (1) 

Py = P2 A UL = uy Xp - us SAY - u 3 
mol) (18) ele (1) (1) (1) 
Ue ee eee See ee ys Qs 
me) 6 Ge) 1) 

ate een ieee 


The same solution expressed relative to the coordinate axes 


OX 41%K+1 15 
= 5 (3.23a) 
> u, {cos6 +e, [C,, +(ik,)(C, Hate Yea toy le 
Oo ih Z 3 
a i u,{-siné te, [Cc +(ik,) (C, Ket ey Yung Tey eld espe Ashes 
e) ne 2 2 
oe p,tl te, {C, +(ik)) (C, X41 *C, eu os ) 1} (3223¢} 
Oo IL 2 3 
Pama ett te, [C_ +(ik,)(c_ x, tC, yj tC, d1} 2 3d) 
O It 2 3 
with 
_ (0) = (1) 
oe = cos 843 ; oa = cos 9 U3 P 
O 3 
a Cl) (ak) ere Ck) 
4 = cos 9, (cos 844, + sin 8 U5 sin 8 3V5 ) 
Zz ee (1) (cls) CL) 
zo, = cos 84! sin 844, + cos 8 U5 + sin 9 V5 ) 
- (0) a exe Cr) 
eS = sin 9 43 P Cc. sin 8 443 
Oo 3 
a 7 1) re (CE) (1) 
cy, = eatsyl a4l cos 844, sin 8 4U5 + cos 845 ] 
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eee ee (1) 2 (1) 
7, Sin 8 4U, sin 8 4Cos 8 U5 + cos 93V5 
- (0) = (1) 
Oo 5 
, (1) : . (1) 
“p, YM, cos 9 4Py ’ “Ds he sin @ 5p ; 
(0) (1) 
C = Mr F C = Mer 
ry k 3 63 Ke 3 
a: (1) ae (1) 
Sa M, (cos 8 ary siné.r, Vine 
= (1) (1) 
=, = M, (cos 9 4X5 + sin 94ry ae 
84 = (F577) /n 


We summarize the procedure that should be followed 
to find the flow field quantities in Region BL, =D below. 

(a) Assume that the flow field quantities in Region B have, 
relative to axes O'x'y' attached to body I, the form given 
by relations (3.3) and (3.4). 

(b) Express these quantities relative to axes OXY) as 
described in Section IV.A.1, by relations of the form given 
by equations (3.6), with 94 = (¥5-¥))/n. 

(c) Use these quantities as upstream quantities for the 
first problem considered and find the perturbed quantities 
in Region Bie as described in Section IV.A.4. The solution 
has the form given by equations (3.22) with k = 1 and the 
reference steady flow quantities involved are given by equa- 


mons (3.17) with k = l. 
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(d) Express the solution relative to coordinate system 
OX5Y5, by using equations (3.23) with o 5 ~ (Po-¥,)/n. 

(e) Repeat step (c) to find perturbed quantities in Region 
B. and repeat step (d) to express them relative to coordinate 
axes OX3Y3- 

(f£) Repeat step (e) till the flow quantities in Region 
BL are found. These quantities are the required quantities 
in Region D. 

The procedure given is well suited for computer 
applications. The program can be set up in a fairly easy 


way using the formulas presented and the number of iterations 


can be increased to the accuracy desired. 


C. AN ALTERNATE APPROACH TO THE EXPANSION SIDE PROBLEM 

In the last section the flow field quantities over the 
expansion side of the wedge were found by a series of itera- 
tive calculations. The same technique can give the flow 
quantities at any point (r,9) in the expansion fan region 
(Figure 3.10). 

In this section another approach for finding the flow 
quantities in the expansion fan 1s suggested. The approach 
consists, as before, of the following steps. 

a) Expressing the governing. equations of motion in the 
expansion fan region in polar coordinates. 

b) Perturbing the equations of motion. 

c) Forming systems of equations that are to be solved 


subject to appropriately formulated boundary conditions. 


147. 





U 


le 1 My, 


a 


Figure 3.10. Illustration of polar coordinate system 
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a 


ds. Governing Equations of Motion 


The governing equations of motion were derived in 


Section II.A.2. They are restated below. 


+ovVv = 0 


4S 
ct]oO 


RE YP = 0 
e) 


abet 


fe 


Er 
e_ - 2 
ie 

au 


O 


In polar coordinates with 


«= =[S=(rA) + mad 

a ae a ee y+ 8 2 
these equations become 
Elbe + Or(V,) +0(V Ja] +O (V,) +0(Vg) 4 +09 (Vg) 


H] 


pir | ie) + (V__) (Val - (Va) [(V_), - (Va) }} 


O{LUV,)o+ (VW) (Vy) al + (Vg) [(VQ) 9 + (VIII 
rent los + (V_) ox] oF (Va) 0g} ea “5 (V_) pol 


x (Va) Pg} ee) 
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(3.24a) 


(3.24b) 


(3224¢) 


(3.24d) 





where p, 0, V. and Ue are the density, pressure and the r- 
and @-components of velocity relative to the polar coordinate 
system chosen (Figure 3.10). Bars over these quantities were 
introduced to indicate their unsteady nature. Quantities t 
and r respresent time and radial distance respectively and 
bars were introduced to indicate that these quantities are 
dimensional quantities. Subscripts t, r and 9 denote partial 
differentiation. 

We assume that the approaching free stream is oscilla- 
Peony. Let . be a small parameter characteristic of the departure 
of the free stream from its mean (average) constant-state 
flow. We also denote by a 0, V.. and Ue the reference steady 
flow quantities (when : = 0) in the Prandtl-Meyer expansion 
fan along a ray 98 = constant. We may express the flow quan- 


tities in the expansion fan as power series of ¢€ in the 


following form 


“aw 


~, 


D = ptept... (3.25a) 
D = ot 25 re (3,255) 
m= V+ evi +... (3.25c) 
Vo = Ve + eV, ere? (3325) 


~e 


where the reference steady flow quantities p, o, ee and Mis 


may be found from relations given in Section II.A.4.b. 


ibs 0) 





We define non-dimensional quantities as follows. 


m= U t/i, i i (3.26) 


oe) 


where 2 is a characteristic length taken equal to one unit 
of length and U. is the average velocity of the approaching 
free stream. 

We substitute (3.25) in the governing equations of 
motion (3.24) and retain only the zeroth and first-order 


terms in ¢€. Equating like order terms in each equation and 


using (3.26) we obtain the following two systems of equations. 


0(Va)_ + Oe(Vp) + 0(V,) = 0 (3. 26a) 
(Ve) L(V) - (V,)] = 0 (3.26b) 
0 (Vo) (V,) + (Va) g] = -P, (3.26c) 
YP Pp = 2 Be (3. 26d) 


r (U9, +0(V_) +0 (V_) 1+ (Vg) g+0 (Vg) +0, (VQ) 


+ p(V_)+0(V_)+0,(V,) = 0 (3.27a) 


Bon (VY) AV (V_) +0 (V,) (CV) .- (V4) Ite (Vv) (VL) 4 


= (V) ItoV, (CV) g- (Vv )] = -rp (3°27) 


8 r 
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rp lU,,(Vg) ,+ (Vg) 2 (V_) 140 (Vg) [(Vg) gt (Vp) | 
+ O(V,) L(V.) +(Vq) 9140 (VQ) L(V) 4(V,) 4] = Pg  (3-27¢) 


r{o[U,p,+p_(V_) I-yp(U,,0, +0, (V_) 1}+(V,) [oBy-YBPq1 


+ (V,) LOD, +OPy-Y (Pe, +pP,) | = 0 (2. 2c) 


Equations (3.26) constitute the governing equations 
of motion in the Prandtl-Meyer expansion fan for steady flow. 
They are satisfied by the Prandtl-Meyer relations given in 
Section II.A.4.b if the additional relation (Va) = a (a 
direct consequemeof the fact that in the expansion fan the 
discontinuities are infinitely weak or Mach waves) is used, 
as can be seen by direct substitution. Equations (3.27) con- 
stitute the perturbation equations for our problem. 


We assume that the unsteady parts of the flow quanti- 


ties are of the following form. 
p> = pympel*t (3.28a) 
9 = oMRet<* (3.28b) 
v. = (¥.)uerXt (3.28¢) 
9 9 
v= (v_)vetkt (3.28d) 
1 bg 


ithe real parts of the complex expressions are considered 


Only. 


2 





where P, R, U and Y are time-independent quantities to be 


found and k is the so-called reduced frequency defined by 


k = o- with w the circular frequency of the oscillatory 


oo 


upstream flow. We also assume that the amplitude and fre- 
quency of oscillations are small and express the time-indepn- 


dent quantities P, R, U and Y as power series in (ik) of 


the form 
p = p92) 4 (ixyp't) 4 (3.29a) 
no Kes (ikyR‘?)? + (32295) 
Swe Se Gayot? y (3.29c) 
ee ec) (3.294) 


We substitute equations (3.29) in equations (3.28) 
and the resulting expressions in the perturbation equations 
(3.27). Retaining terms of zeroth and first-order only and 
equating like order terms in each of the equations we obtain 
the following two systems of equations. 

zeroth order equations 


(0) 
0 


(0) 


(0) 
JR al 


oM[(V_)+(V_), 4. eMI(V,)R tr(V_)R 


. a SACU) eee) eer (0) (0) as 
xi (V4) lo,U +0U, J+o(V_) [¥ aa ] = O (3.30a) 


SS 





o(V,)ME(V_) 9- (Vy) JRO”) +rpyMP t°) 45 (V5) L(V) 2-(V,) 10? 


21, (0) (0) ee 


-0 (V, ) +0 (V_ ) L(V, Yo 12s Ee yx Oe 5.0 OE ) 


O(V,)ML(V_)+(V 9) gIR 9) +pyMP 6°) +0 (V,) L(V) #(Vg) 910 


om (0) ee Oo 
+0(V, ) trp(V_) (V,)UL to (Vi) WV A)Y = 0 (3.30c) 
= (Vg) Dy “Mo ,P °) +opyM I (V,) Po ar (VW pee 
-(¥,) oMp,R'°) -ppyM(V,) Ryo +r (V_) RO} 
+ (V,) [op,-ypp,Ju'?? = 0 (3.304) 
9 9 6 


First order equations 


ae ee a) ns 5 CL) 
OM (V,) H(V_) IR +oM[(V,)Ra° tr (V_)Ro°" | 


+(7,) (00,7) 49.07) 140 (7) ty) +ry 2) | 


= -rpmu_R‘°? (3) 3ia) 


re ee ere) lye) yg 
O(V,)ME(V) .-(V,)IR© +rpyMP ”’ +0 (V,) [(V_), 


(1) Zool) C5) 


(v,)]u‘*) 5 (9) “u'*) +0 (V_) L(V, ) ¥e 
tr(vy¥)) = -rp (7 yuiy (3.31b) 
3 ee aa ere 
o(V,)ME(V_)+(V,) IR" +pyMP,” +0 (V,) L(V) 
as (1) ee ee (te 
+(V,) ,10 +o (V,) a eee) 
+0 (0) (¥,)¥7) = -rpu,(v,)u'°? (3.31c) 
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2. Boundary Conditions 


The proper boundary conditions to be prescribed are 
found by realizing that the problem is an initial value problem, 
1.e., the upstream conditions completely determine the solu- 
tion and the flow matching conditions must be imposed at a posi- 
tion which differs by O(« from that of OB and hence must be 
determined as part of the solution. 

oe oOlution 

Because of time limitations and unsuccessful choice of 
test solutions we have not been able to find a solution. 

We should note that once a solution is found (either 
in closed form or by use of a computer program) the flow quan- 
tities on the upper surface of the wedge (Region D) can be 
obtained by using the method described in Section III.A.2. The 
same method may be employed to extend the solution of the prob- 
lem to the case of an oscillating wedge in oscillatory flow, 
since in the last case, the oscillations of the body are not 
expected to influence the flow field in the expansion fan 


region, except in the neighborhood of line OC. 


ILS 





V. SUMMARY AND RECOMMENDATIONS FOR FUTURE WORK 


A. SUMMARY 

The perturbation method proposed by Professor Hui, as 
applied to a two-dimensional, oscillating flat plate, with 
attached shock wave, at an arbitrary angle of attack, ina 
steady, inviscid supersonic/hypersonic flow, has been des- 
cribed in detail. For periodic oscillations of small ampli- 
tude and frequency, the first-order closed form solutions for 
the flow field quantities in the disturbed regions, over 
both sides of the two-dimensional flat plate, have been given. 
The “in-pitch" stability derivatives of the oscillating flat 
plate can then be predicted and a criterion for the neutral 
damping boundary can be obtained. The flat plate results 
may be naturally extended to include bodies of slightly more 
complicated shapes and this was done by Hui for Caret wings 
and wedges [{Ref. 7]. Utilizing strip theory these results 
may also be extended to three-dimensional bodies of a similar 
cross section and this case was described in Section IIl for 
a flat three-dimensional wing of arbitrary planform [Ref. 5]. 
The “in-pitch" stability was studied and the agreement with 
potential flow theory results (for zero angle of attack) 
was found to be good. 

The extension of the above results to bodies with cross 
sections composed of curved segments or straight segments 


forming downstream corners, is possible (Figures 4.1 a,b). 
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In Section IV the perturbation method proposed by Hui 
was extended to include the effects of periodically oscillating 
upstream flow over both sides of a stationary flat plate. 
The method may be readily extended to include oscillations 
of the plate. The assumed form of upstream flow oscillations 
includes as special cases the form of induced oscillations 
in the flow field by the oscillating flat plate considered 
before in steady supersonic/hypersonic flow. The solution 
over the compression side is given in closed form, while 
that over the expansion side is given as a series of iterative 


calculations. 


B. RECOMMENDATIONS FOR FUTURE WORK 

The upstream oscillatory flow solutions may be combined 
with the upstream steady flow solutions to study three 
general problems. 

First, the induced flow field by two-dimensional bodies, 
of any cross section, at arbitrary angles of attack and 
attached shock waves, in steady flow. For this problem the 
flow field can be obtained in steps, from leading to trailing 
edge, over the straight line cross section segments (Figure 
4.la). Curved segments could also be approximated by a 
series of straight line segments (Figure 4.1lb). 

Second, the induced flow field by a formation of bodies 
with the exception of 

i) Regions behind crossing shocks, 
ii) Expansion fan overlapping regions, 


iii) Mixed regions of above. 


Sy, 


= 








(a) see 


REGIONS 8B 





(b) 


Figure 4.1. Illustration of flow solutions 


Jepeyts. 





Third, the induced flow field behind three-dimensional 
bodies of any cross section, utilizing the strip theory 
approximation. 

The dynamic stability, in pitch, of the bodies involved 
in each of the above general problems, might be studied. 

A closed form solution for the expansion-fan problem 
formulated in Section IV.C might be sought. Similarly 
a closed form solution might be sought for the more general 
expansion side problem considered in Section IV.B. This 
solution might possibly be found by letting the step angle 
8, (representing the step change in flow direction) tend 


d 


to zero in the basic formulas that are iteratively used. 


og 





APPENDIX A 


PROCEDURE FOR SOLVING SHOCK BOUNDARY CONDITIONS 


Let 2, = Dh careele 
X= wee, Vy - gradG} 
at 
X. = pee 7, * gradc] 
gt 


where & is the characteristic length chosen, G = 0 is the 
equation of the shock and the subscripts 1 and 2 denote 


quantities before and after the shock respectively. 


Equations (2.3a-2.3c) are written with u = a 
PX, = P4X5 (A. 1) 
2 2 2 2 
0,X) + P,z = 5X5 + poz (Ax 2) 
p ss 
x? + u + z° = x? + u 2 5° (A.3) 
1 Py Z P5 


Consider a harmonic oscillation with frequency w and maximum 


angular deviation of the unsteady flow from the steady reference 


Flow e. 
femec = Panes where 

tee 
t (non-dimensional time) = —t 

X 

wd 

k (reduced frequency) =e 

fo) 
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t: dimensional time 


uo: reference velocity (in steady flow) 


For the general case assume that the approaching flow is also 


unsteady and let 


x = Xia = EX, 5 Xo = Koy + EX, 5 zZ = 25, + &2,5 
fees Fir. ~ *P12 Po Sa ey, 
ge, 681 «6K P02 en oneg 25 


We substitute these values in (A.1-A.3) and retain only zeroth 
and first order terms in ¢. Equating zeroth and first order 


terms we get the following two systems of equations respectively. 


P54 %91 = 057401 (eee) 
Po, [e54,7;, (lw) = yy fey, +92, (1-H) J Uevo Beal) 
o> Lp x? tip ..2.,}) = to) as x? (ae a) 
eg a1 HP 711 ae ae 
aa oe ein io 7 Palo“ © 9 20“91 (A.1.2) 
Poo (P547KP qq) = Py 2 (Py tkP oy) —Po0 (Py ktP 91) tPy 9 (P44 ~KP 9) 
; (A.2.2) 
Poo (Py. XqythPy 7247) = 204% yy Xp9l (url! 4-97) 
2 
ee le een "on Pi 212°P 12711: a! 


with k = u-l1 = (y+1)/(y-1). 


Ga 





Equations (A.1.1-A.3.1) give the unknown flow quantities 
behind the discontinuity, for the steady reference flow (in 
terms of known flow quantities before the discontinuity) and 
are also used to simplify equations (A.1.2-A.3.2) even further. 
Pieseequantities will then be used in simplified equations 
(A.1.2-A.3.2) to give the unknown flow quantities behind the 
discontinuity for the unsteady flow. 

To complete the system of boundary conditions across the 
shock we now consider equation (2.3d) which expresses the 
conservation of tangential momentum. This equation is 
equivalent to 

ee 2 =) VO Tt (A.4) 
where T is a tangent to the surface vector, such that 
VG°T = 0. 
For a surface G = xtango ~y +€Q(x) = 0 we will choose 


{l,tan @+¢eQ'}. 


this vector to be T 
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APPENDIX B 


OSCILLATING WEDGE IN STEADY UNIFORM FLOW 


For this case which was considered in Section III.A.l 


Pumeuce 2.la), 


< 


i= ae Usicos ¢)7-sin 6} 


<| 
il 
at 
& | 
<| 
a] 
II 


uo il ety. eV } 


Z Zee 
G = x tango -~-y +€Q(x) 
Then 
7G. = Lftang EO path De =e {1,tan 9 +eQ,} 
ey 
Z° = @ + €2Z = Jl + Sem 6 + e[2 tan?oQ_] 
11 12 x 
faee= X_.. + &X = (2lUeeee (cos eo) 4 Oi K)Q]} 
1 ll eZ e°cos ¢ Xx We 
xX, = Ko, + €X55 = uo{tan ¢ + e(Utando -V Beet (ik) Qj 1 
We also have 
ees Py) we Pathe OPES Poy TEP 22 = Po TEPQIMS® 
Pee «Pa +t FPyp = Po FEO 95 = 091 EP og *PQ FEPQLMLR 


G3 





These values will be substituted in equations (A.1-A.3). 


Equations (A.1.1-A.3.1) are considered first. From Equation 





Z 
(A.3.1) we get after multiplying it by gos 
OU. 
p Y Pp 
S| ee 2 co Vad ere 
S—(Sa5N es cs Sa } = -—~< sins (eerie ls) 
O ie ous y-1 
: Y Po 1 P56 ; 
Setting — > and solving for a get 
26 (y+1)M2sin*3 
, a ao (Aron b AL) 
co (y-1)sin BM, +2 


See@eequation (A.2.1) we get after dividing it by op, and 


. Po 
Gemving for 5. 


oo 


O 
1S, 


Po y+1 O y+1 


10) 

8 
© 

8 
a2 
fe 


Substituting for 2 67? oo its value from (A.3.1.Al), we get 


a 2ysin“gM* - (y-1) 
= = a (Av2elsAl) 


From equation (A.1.1) we get after multiplying it by cos 6 
euecriee = O44 sin db CAP eve) 


From equation (A.4) we get the following two equations equating 
the zeroth and first order terms in e« and using the relation 


cos @ cos ¢ - sin @ Sin 6 = cos 8 


UW eS (b= u, cos ) (A.4.1) 
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-U,sin 6Q, - u, (U + tan $V) (A4% 2) 


From equations (Al.1.0) and (A4.1) we have 


P65 _ tan 8 


Q tan ¢ 





(Feist. Al) 


Mmemations {(A.1.1.Al), (A.2.1.Al) and (A.3.1.Al) consti- 
tute the supersonic uniform wedge flow solution presented in 
section II.A.4.a. 

Equations (A.1.2-A.3.2) are considered next. The time - 
independent unknown functions R, P, U and V can be found in 
terms of Q' and (ik)Q in the following order. From simplified 
equation (A.3.2) the function R; from the known function R 
and equation (A.2.2) the function P; finally from equations 
(A.1.2) and (A.4.2) the functions U and V. The derivations 
are lengthy and tedious and will not be included here. The ex- 
pressions for the complete set of coefficients in equations 


(2.11) are given in [Ref. 7] and are quoted below. 


~~ 


Q Q Q 
Meee eos ¢ [142 (M-=1) tan-¢-yw(—> -1)] (Al. 3a) 
06 2, © 0 6 
ui P 0 2 Po Po 
(l-w)B = (1-—)cos"o[1l +— W-yWw(—-1)] Nes Sie) 
3E a a 
a P 4 y-1.,,°0 
(1l-W)C = 2kH(1 -<—) cos o{l To ee (Al. 4a) 
O fore) 
D = Oy eae (A1.4b) 
E = (A-B) eon 0 (ar Say) 
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F = --B tan 6 (Al. ob) 


E ee 4 Po 
(l-W)G = Ta) ON Giese el (Al .6a) 
O oO 
: ie Z 
J = G/cos¢ (Al .6b) 
M 
Kk = aoe ee H = Maeeol: tan ¢, W = we in Co 
omar fe) Oo 


fo) 

To illustrate the derivation procedure suggested, we 
find, below, the expression for the function R and compare the 
coefficients G and 3 with those given by relations (Al.6a) 
and (Al.6b) above. 

We multiply equation (A.3.2) by cos*$/p_U and use relations 
Memes. Al), (Al.2), (A.4.1), (A.1.1.Al) and cos 9 = cos (8-6) 
to simplify it. We finally have 
2 cos (pe Wel: 


ee 2 
M, sin $ 5 i: eae: + (1k)Q/cos 96] 


~~ Stead 


Thus the coefficients G and J found are 


os 3 0 
2 cos’ 6 00 yol 

GC =F a (— - SO) (Al .7a) 
M, Sin 60 a ytl 

~s ~~? 2 

J = G/cos”6 (Al .7b) 


~~? 


The second coefficient J is expressed by the identical rela- 
tions (Al.6b) and (Al.7b). It can also be seen that the 
expressions (Al.6a) and (Al.7a), for the first coefficient 


G, are identical by direct substitution and use of the 
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following additional relations given in Section II.A.4.a, 


Po = 2 + (y-1)x' won 2 eis 
Po yay" : ZX” = Ge) 
where 
a. — M“sin“8, 
x = m*sin*¢ 
Oo 
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APPENDIX C 
OSCILLATING FLAT PLATE IN STEADY UNIFORM FLOW-- 
EXPANSION SIDE 
For this case, which was considered in Section III.A.2 
(Figures 2.2a and 2.2b) we assume that there is no change 


in the steady flow quantities across the front. Thus 


V5 = u, {1,0} ; Us = u, {1 +eU,ev} 
G = xtang -y +€Q(x) 
Then 
VG = titan do + €Q',-1}, Tee cand. + eo! | 
ie 
eee ez = 1+ tan*¢ + e[2 tanoQ'] 
i 12 
= — | " 
X, = X,, + &X. = u, {tan 6 + e(Q' + (ik)Q]} 
X, = Xo, + EX, = u, {tan d+ e(U tano- V + Q' + (ik)Q]} 


We also have 


a we a EPin = Py + €°0 
Dy = Dag 7 ED ee mesh Sh eens 
By Ba oP Sey eS 
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+ ep, mM R 


25) i eee Py il 


Equation (A.3.2) is simplified using equation (A.3.1) and 


becomes with Pio = OF 


e) 
s > ied a 
P52 = Po, (1 y+ eZ 


Subsituting values given above we get 


Ree GO’ vr J (1k) 0 (AZ) 
with 
4 see" Z 
G = (y+1)M, sing ’ J = G'/cos ) (CAP eel.) 


Equation (A.2.2) is simplified using equation (A.2.1) 
and beomces with P49 = Pi2 = 0 


Pa Saabs ee! 


0 p.. (0..K-0..) 
22 Py (P51 K — P44) 


P52 


We substitute values from above and have 


P = R (AZ. 2) 


From equations (A.1.2) and (A.4) we get the following two 


relations respectively, 


U tan o@ =- V = =M, tang R 


U+Vtan 6 = O 
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Solving these relations we get 


V ins A'Q' + B'(1kK)Q tAZ, 3) 
with 
a = 4cos"o At 2 AD fase eee) 
Wai 1 = cos o 35 
and also 
UPe= EE Ole F' ak) oO (A234) 
with 


Nestor sin 6 2 
E' = = ears aia : F' = E'/cos 6 (A2.4.1) 


Equations (A2.1-A2.4) give the boundary conditions 
weebor=2.161) and (2.17£-2.171) when the power series expan- 
Sions of the time-independent quantities, given by relations 
(2.8), are used. The apparent differences between the ex- 
pressions given here for the boundary conditions and those 
given in [Ref. 8], formulas (21) and (22), are due to the 


different representations of the equation of the front. 
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APPENDIX D 


OSCILLATORY UPSTREAM FLOW PAST A STATIONARY WEDGE 


For this case which was considered in Section IV.A.3 


(Figure 3.4b) we obtain from relations (3.6) and (3.9) with 


E5 = 0, 
Vy = a 05 7&1 ey eee 6, +€,Vz} 
Men ie P12) Po.” “1Potes co 
Mm - "11 78y2 = %% +6105 Rp (A3.1b) 
V5 = uf) +¢,0,,€,V,} 
pe 2 P22 = Pe Polo 1 aes! 
eee ee “1929. ~ 26 =i? ot ae 
where Up, Vas P, and R, have the following general form given 


by (3.6a-3.6d) 


M, = Ch + (ik,) (Cc, x+C_ y+c. ) CA30)) 

O i: 2 3 
We also have with La = 2 and G = x tang -y +€,Q, (x) = 0 
ae. zl ' o nm 1 

VG = tan ¢+€,Q0), Lee co slit b+ €,Q7} 
Q 

ie = @ +e.Z = L eet Me -2 tanoQ, (A3.1le) 
11 ele If Ht 


bby 





a = i" B/eos ¢ (A3.1f) 
Xo = Yo. Up tan $-V, + cos 2 aoe Meat ili (A3.1g) 
Koy = uy tan? (AS 21) 
Ko5 = ujlU, tan¢ ~V, +Q) + (1k) )Q)] (A311) 


Equation (A.3.2) is simplified using equation (A.3.1) 


and becomes 


2 
1 ®an, Sipe Soon 
7 jo.) — 

11 lak P31%11 


Using relations (A3.1) we obtain 


v9 
ll 


; 
l KR OB *KR Ye rag + RRB +R 9] 43 (ik)Q) (Aeia2a) 


with 


aR, = Hy Sind/sin 686, *R, = -Hy cos ¢/sin 8, 
Kp = 1/M +H, /2, Ke = - —____, +, 
3 . 4 M [sin“8(y-1)Mo +2] 
© O 
B 
a 2 = 
BR. = H, cos O7 tare cee PR = H, /tan B 
e) 
Bd. BOs 
al 2(1 y+l 0 ) M4 
Op 


IL 





Equation (A.2.2) is simplified using equation (A.2.1) 


and becomes 


ee 2 Pi Pa! Pan {Pit ten}) Pa! 
D.. [po..- (unL)o..) Da ieee—=tu-lNons 1 
p 
De 
ae 
aie (pa 


We substitute values from (A3.1) and get 


P = K U_+K V_+K R._+K P_+K '+kK ak A3.2b 
p> C P, B Py 8 pee Be 1) 24 ( ) 


mR = Hae, 45 1,9 .856 
1+ (u=D)P,/Pg_ 

K = HK + 

p;  "2*R, * “MOTE -W=1)0,70, T 


eS 


vi = ‘x, + 1/ (YM). 


Ah Sires eae) 
eee eo Tl =a oe oe 


SB B 


From equations (A.1.2) and (A.4) we get a system of two equa- 


tions and solving for Uy and Vy we obtain 


Vv, = Ba hy Tey Paty Ory 8 (aa2c) 


Lye 





vy SOR 
4 J 
K = HK + H,, K = HK Hoy tan 6 
v5 3 Ry 4 V5 3 R. 4 d 
K. = HK, + Hy, cos(8+o)/tan 6 
5 5 
K = HK, - cos*¢lp_ /p_- 
Ve 3 Re On, 0 
H. = MS tan®o, Hy = sin 8 Sind cos ¢/sin 8 cos 8 
U = allt V..+K +K,, P_+K,, Q!+K,, Q., (1k) (A3 .2d) 
u u u 
1 Uy BU, Be UBT AU, B poets 1 
with 
= = tan oK (j = 3,4,6) 
J J 
K = K tanod+u_ /u_, K = (-K tanod+u_ /u_)tan 6 
Uy Vy oO, 0 u, V5 On 0 
K = —-(K tang@d+u  siné,/u_) 
U, Ve Op a =o 


The general form of the shock boundary conditions expressed 
by equations (A3.2) is 

— + Bt ° 
4 K oes VerK KY TVS O! LG, Q, (1k, ) (A3.3) 


i vy Yo B Yee 


Using relation (3.10) the zeroth and first order boundary 


conditions across the shock become 


0) (0) (0) (0), ae 
Yy = Ky Up *Ky Ye oe * fae gi (A3.4a) 
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J (il) (1) (1) (eis) ol) (0) 
Yi = Pee +K. V *Ky Rp ve SI +K,, QO 


Y) B Yo B B 6 iL 
(A3.4b) 
Introducing (A3.0) we obtain 
my” = C and ma? = C x+C eC 
Mo we 2 3 
Relation (A3.4a) then gives 
co) _ , (0) Bs 
Yy = se igs + Ky At y = xtan 90 (A3.5a) 
where 
K <0 ee ae se + Gs: + K C 
a ino to — Xo Eg Po 
K = 
yo" %, 
Similarly, relation (A3.4b) gives 
(1) an ia) 
Yy = K. x +Ky, y+K, +k, Q) +K, Q, (A3.5b) 
XX VY x Oo 
At y = Xx tan 9 
where: 
K = K., C + KC + c ~ C 
Y ex Sey al se SL we oT 
K = K_C ares (0: + & thee 
Toy vy us Yo V5 Me ro vy D5 
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APPENDIX E 


ENSPITCH DERIVATIVES BASED ON TN*S 2699 AND 3196 


To find the stiffness and damping-in-pitch derivatives 
given in [Ref. 5 ] we proceed as follows. 


We define the following quantities: 


S = planform area 
b = half wing span 
Cc. = aafelere (elelonatsl 
Cc, = jspijs) (es ele pao 
Lie = leading edge angle of sweep 
M = free stream Mach number 
B = cotangent of Mach angel = ey oe 
A = aspect ratio = (2b) */s 
XK = taper ratio = c,/C, 
C = mean aerodynamic chord = 2C_(A*+A41) /3(A+1) 
m = slope of LE = cot Lie 
A' = parameter = AB 
m' = parameter = mé 


We also state the following stability derivatives defined 
eeqolicitly in [Ref. 11,12]. 


Resulting from steady state motion 


L = lift curve-slope derivative 


CS = lift and pitching moment derivatives (due to 
gq Mg steady pitching) 
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Resulting from time-dependent motion, 


lift derivative 


Q 
ul 


Q 
il 


pitching moment derivative (due to constant 
Q vertical acceleration) 


To account for differences in definitions of quantities 
we form the following factors 


el _ Z 


The stiffness and damping-in-pitch derivatives as defined 


in this thesis are then given by 


& = FC ’ Core a oe eet Ce) (ejse) 
9 a m9 2 ke mo 


To find Cn ; Ca , ox from NACA TN's we proceed as follows: 
o q O 


meer orm 6, A, A’, m', ee al: 


ii) Read from NACA TN 2699 (Figures 11-15) BC and find 
oi 
Ci 7 ((8C, )/8 
oi a 
iii) Read from NACA TN 2699 (Figures 26-30 8C and use 


L 
Q, 


its value to find pivot position as follows 


ACW a? een) (cea) 
fo Q 


iv) Read from NACA TN 2699 (Figures 16-20 and Figures 21-25) 


BCE and 8Co which refer to pivot position found above. 
q q 
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v) Use following formulas to find BCL and BC. 
gq q 


‘@) 
I 


(BCE - 28C_ )/8 
q gq O 


Q 
il 


[scx + (BC ) (BC, )/(8C, )I/B 
q S| 0 q o 


vil) Read from NACA TN 3196 (Figures 11-15) quantities 


(BC, .), and (8C,.) » 


a a 
follows, 


and use them to find Cre as 
Ot 


_ 2 
ci. = (M*(8C, 


& 


2 eu ye 
) 78" + (HT7B" +2) (BC, .) 91/8 


vil) Use relations Bl to find comparison values for Ch, 


and eas 
6 


di ie, 


G 





JL)e 


a. 


IL eae 
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